PERTURBATIONS OF ORTHOGONAL POLYNOMIALS WITH 
PERIODIC RECURSION COEFFICIENTS 



DAVID DAMANIRi, ROWAN KILLIP^, AND BARRY SIMON^ 



Abstract. We extend the results of Denisov-Rakhmanov, Szego-Shohat- 
Nevai, and Killip— Simon from asymptotically constant orthogonal polynomials 
on the real line (OPRL) and unit circle (OPUC) to asymptotically periodic 
OPRL and OPUC. The key tool is a characterization of the isospectral torus 
that is well adapted to the study of perturbations. 



1. Introduction 

This is a paper about the spectral theory of orthogonal polynomials on the real 
line (OPRL) and orthogonal polynomials on the unit circle (OPUC), that is, the 
connection of the underlying (spectral) measure and the recursion coefficients. 

Specifically, given a probability measure, dr], on R with bounded but infinite 
support, the orthonormal polynomials, Pn{x), obey a recursion relation 

xpn{x) = an+iPn+i{x) + fe„+ip„(.T) + a„p„_i(a;) (1.1) 

where the Jacobi parameters {a„,6„}^i obey bj e R, aj G (0,oo). As is well 
known (see, e.g., [1021 Section 1.3]), (|l.ip sets up a one-one correspondence between 
uniformly bounded {an,fen}5^=i a-nd such measures, drj (this is sometimes called 
Favard's theorem). 

Similarly, probability measures, dfi, on 9D which are nontrivial (i.e., their sup- 
port is not a finite set of points) are in one-one correspondence with sequences 
{a„}^o of Verblunsky coefficients in D = {z : \z\ < 1} via the recursion relation 
of the orthonormal polynomials (y9„(z), namely, 

Zipn{z) = pn(pn+l{z) + a„(/?* (z) (1.2) 

where 

<p*Jz) = z-^^iT/^ P„ = (l-|a„P)^/' (1.3) 

Underlying the association of measures and recursion coefficients are matrix 
representations. For OPRL, we take the matrix for multiplication by x in the 
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{Pn}^=o basis of L^(M, d?7), which is the tridiagonal Jacobi matrix 





bi ai 





J = 


ai 62 






02 






V 





(1.4) 
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For OPUC, one takes the matrix, C, for multiphcation by z in the basis obtained 



by orthonormahzing {1, z, z ^, 



} in L^{dB),dii). This CMV matrix (see 



[I02I Section 4.2]) has the form 

C = £M 

C ^ eiaa) (S e{a2) ® ■ ■ ■ 
M ^ lixi©e(ai)®e(a3) 



eia) 



(1.5) 
(1.6) 
(1.7) 

(1.8) 



where p = (1 — |ap)^/^. Note that C is unitary, while J is self-adjoint. 

As a model for what we wish to prove, let us briefly survey some of the main 
results relating to (slowly decaying) perturbations of the free case, that is, a„ = 1, 
&„ = for OPRL and a„ = for OPUC. 

(1) Weyl's Theorem [Hll O [H |45l H] . If a„ ^ 1, &„ ^ 0, then a^ssidv) = [-2,2] 
and if a„ 0, then aessid/J.) = SB. Here aessidij) is the (topological) support of 
the measure, drj, with all isolated points removed. 

(2) Denisov-Rakhmanov Theorem [89l [90l [771 Ell US]- If Ccss(d?7) = Eac('i?7) = 
[—2,2], then a„ — > 1 and 6„ 0. If (Toss(d/i) = Eac(<i/i) = 9D, then a„ ^ 0. 
Here I]ac(rf'7) is defined as follows: let drj — drj^c + drjs with ^773 singular and 
dria,c = fix)dx, then J^s^ddr]) = {x : f{x) ^ 0} as a measure class, that is, modulo 
sets of Lebesgue measure zero. 

(3) Szego's Theorem [lOS [Iini [Ml [H [611 . In the OPUC case, define 



z{d^l) = - / log 

Then Z(dfi) < 00 if and only if 



de 



de 
2^ 



In the OPRL case, define 



Z{dv) 



dE 



Then, if we assume supp{dr]) C [—2, 2], we have 



N 



Z(dp) < 00 



lim sup log(aj 



(1.9) 



(1.10) 



(1.11) 



(1.12) 
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and if that holds, then 



and 



^(a„-l)2 + &2 (1.13) 



^(a„-l) and (1.14) 



are conditionally convergent to finite numbers. 
(4) Killip-Simon Theorem ^61/. For OPRL, define 



Qidv) 



and let {Ej} be the point masses of drj (eigenvalues of J) outside [—2,2]. Then 
(fTTS)) holds if and only if a ess {dr]) = [-2,2], Q{dr]) < oo, and (l^^j I - 2)^/2 < oo. 

(5) Nevai's Conjecture [SI EJ. For OPRL, if J2n=i\^n - 1| + \b„\ < oo, then 
Z{dp) <oo {Z given by (fTTT]) ). 

The five results listed above capture different aspects of the philosophy that the 
measure is close to the free case if and only if the coefficients are asymptotic to the 
free ones. In this paper, we study extensions of all these results to perturbations of 
a periodic sequence of Jacobi or Verblunsky coefficients, that is, 

= = n>l (1.16) 

or 

«i+p = «l°^ n>Q (1.17) 
and some fixed p> 1. Note that p = 1 is the perturbation of the free case considered 
above. For simplicity in the OPUC case, we will normally suppose p is even — indeed, 
the shape of a CMV matrix repeats itself only after shifting by two rows/columns. 
As explained in Section [T51 the situation when p is odd can be reduced to this using 
sieving. For OPRL, p is arbitrary. 

The philosophy described above becomes more subtle when we move to the 
periodic setting; rather than having a single 'free operator' we have a manifold 
of them (the isospectral torus). Nevertheless — and this is the main thrust of the 
paper — spectral measures that are close to those of the isospectral torus correspond 
to coefficients that approach the isospectral torus. One of the key obstructions 
here is that a sequence of coefficients may approach the isospectral torus without 
converging to any particular point therein. 

In order to make these heuristics precise, we need to make a few definitions. To 
keep the presentation as coherent as possible, we will focus our attention on the 
OPRL/ Jacobi case for the remainder of the introduction. 

To any pair of p-periodic sequences, {an\b^''}nez, we can associate a two- 
sided Jacobi matrix Jq. Two such pairs of sequences are termed isospectral if the 
corresponding Jacobi matrices have the same spectrum. We write Tj^ for the set 
of p-periodic sequences that are isospectral to Jq. Topologically, this is a torus as 
explained in Subsection 12. 141 below. 

Given two bounded sequences {a„,6„}^]^ and {a^i,6^}^i, we define 

oo 

d™((a, 6), (a', 6')) = ^ e-'=[|o™+fc - a'„+J + \b^+k - 6;„+fe|] (1.18) 

k=0 



4 D. DAMANIK, R. KILLIP, AND B. SIMON 

which is a metric for the product topology on X^((0, i?] x [—R,R]). The OPUC 
analog is 

oo 

d™((a), (a')) = ^ e-'=|a™+fe - (1.19) 

a metric for X ,^ D. The distance from a point to a set is defined in the usual way: 

dm{{a, b), T) = inf{d„((a, 6), (a', b')) : {a',b') e T} (1.20) 

and similarly in the OPUC case. 

We begin with the periodic analog of Weyl's Theorem. 

Theorem 1.1 (Last-Simon 69j). Let Jg be a two-sided periodic Jacobi matrix and 
J a one-sided Jacobi matrix with Jacobi parameters {am,bm}m=i- V 

d„,((a,6),TjJ^0 (1.21) 

thei^ 

CToss(J) =^t(Jo) (1.22) 

As indicated, this result first appeared in (6^. It is derived from a theorem that 
had earlier been proven with different methods by others [121 EH [EH] ■ The inclusion 
CossiJ) 3 cr(Jo) follows easily using trial vectors; the reverse seems to be more 
sophisticated. In Section [SJ we prove this using the methods of this paper. The 
OPUC version appears here as Theorem 115.11 it was also proved in [53] . 

Note that (|1.2ip does not imply that there is a sequence {(a^,6^)} € Tj^ such 
that 

d„((a,6),(a',5')) ^0 
It is much weaker. Equality of essential spectra under this stronger hypothesis 
follows immediately from Weyl's original theorem on compact perturbations. 
Our first major new result is an analog of the Denisov-Rakhmanov Theorem. 

Theorem 1.2. Let Jo be a two-sided periodic Jacobi matrix and J a one-sided 
Jacobi matrix with Jacobi parameters {am7&m}m=i- V '^cssiJ) = c(Jo) and 

Sac(J) = a(Jo) (1.23) 

then dmiia,b),Xjo) ^ 0. 

Remark. Using Theorem ll.4l below. we will show that the hypotheses of this theorem 
can hold while (a, b) only approaches Tj^ without actually having a limit. 

A two-sided p-periodic Jacobi matrix is said to have all gaps open if the spectrum 
has exactly p connected components — the largest number possible. As explained 
in Section [21 this holds generically (indeed, on a dense open set). 

Our next new result is 

Theorem 1.3. Let Jo be a two-sided periodic Jacobi matrix with all gaps open 
and parameters {an\bn^}'^^_^. Also, let J be a one-sided Jacobi matrix with 
parameters {an,6n}^i and spectral measure drj. We assume that acss{J) = cr(Jo) 
and 

J2 dist(^j, (Te.s( J))'/' < OO (1.24) 
j 



Recall that acss(J) is obtained from the spectrum of the Jacobi matrix J by removing all 
isolated points. 
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where {Ej} enumerates the eigenvalues of J outside cr(Jo). Then 



log 



(^) ^ ^ '^(■^o))"'^' dx<oo (1.25) 



implies 

pN . . 

exists and lies in (— cxd,cxd). Conversely, (jl.25p /loWs so Zoni? as 

limsup Vlogl I > -oo (1.27) 

and in this case, the limit in (jl.26p exists and lies in (— cx),cx)). 
Lastly, if (fT^ or (fTTT)) /loZrfs, i/ien 

OO 

Y,dm{{a,b)Xhf <^ (1-28) 
anrf there exists Ji G 7}^, so i/iai 

d™(J,Ji)^0 (1.29) 

Remarks. 1. Thus, when (|1.24p - p.27p hold, J has a limit Ji in Tj^. In the normal 
direction to Tjg, the convergence is £^ (in the sense of p.28p ). But in the tangential 
direction, we only prove it has a limit. It would be interesting to know what can 
be said about how slowly p.29p can occur and to know if there are examples where 
(fmD ^ lfr^ hold but 

OO 

Y,d^{J.Jif = ^ (1-30) 

Tn=0 

2. Notice that (|1.27p will only fail if the partial sums converge to — oo. 

3. The final statement that there exists Ji G Tjf, with (|1.29p is not our result 
but one of Peherstorfer-Yuditskii 'Wf\. With our methods, we can prove that the 
a's and 6's approach a periodic limit only if we replace (jl.24p with the stronger 
assumption that the discrete spectrum is finite. 

4. By (11231), aU af^ in fL^ and PTTTI) can be replaced by Cap(o-(Jo)), the 
logarithmic capacity of the spectnmr of Jq. 

Our third new result is 

Theorem 1.4. Let Jq be a two-sided periodic Jacobi matrix with all gaps open and 
parameters {an\bn'^}nei.- Let J be a Jacobi matrix with parameters {a„,6„}^]^ 
and spectral measure drj. Then 

OO 

^ d„((a,6),rjj2 <oo (1.31) 

m=0 

if and only if 
(i) acss{J) = a-(Jo), 
(u) ^ dist(i?j, CToss (-'")) < OO, and 
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(iii) - / \og(^]dist{x,R\<j{Jn)y^^dx <oo. 

Ja(J„) \ dx J 

Here {Ej} enumerates the (discrete) spectrum of J outside cr(Jo). 

Remarks. 1. Since (i)-(iii) are equivalent to (|1.3ip . one may easily construct exam- 
ples where (i)-(iii) hold, but there is no Ji with (|1.29p . This provides the examples 
promised in the remark after Theorem 11.21 It also shows a stark difference be- 
tween (|1.24p - (|1.25p and (ii)-(iii). In terms of the spectral measure, this difference 
is reflected only in the behavior near the band edges. 

2. As we will see (Section [T^ . there are results even if all gaps are not open, but 
for Theorem 11.41 they are not so easy to express directly in terms of the a's and b's. 

3. A special case of part of Theorem 11.41 is known, namely, Killip [53] proved 

that Sac(J) = cr(Jo) for {a„,6„}^^o obeying ^^Ja„ - al°^|2 + |6„ - < oo 

(which is a strictly stronger hypothesis than (|1.3ip ). 

Theorem 1.5. Let Jq be a two-sided periodic Jacobi matrix and J a Jacobi matrix 
with Jacobi parameters {an,&n}^i md spectral measure drj. Suppose 

oo 

^|a„-a(°)| + |6„-6(°)|<«3 (1.32) 

n=l 

Then (fT:25)) holds. 

Remarks. 1. Condition (|1.32p can be replaced by 

oo 

^d„((a,6),r7j < oo (1.33) 

n=l 

Indeed, if (fT35|l holds, then holds with {al^\b';°^ } replaced by some fixed 

sequence in Tj^ . 

2. As we will show (see Proposition 13. Sp . the theorems above continue to hold if 
dm is replaced by 

p-i 

dm{{a,b), {a',b')) = ^(|a„+fe - a'^^^] + |6m+fc - 6™+fc|) (1-34) 

fe=0 

For OPUC, we need to sum k from to p in order to get an equivalence; see the 
discussion at the end of Section [H 

In Section [151 we prove an OPUC analog of each of these theorems. We need to 
replace the all-gaps-open hypothesis with a stronger one (that holds generically) . 
The deficiency is not so much with our method, but rather that an independent 
question, which is known in the Jacobi case (independence of the Toda Hamil- 
tonians), is currently unresolved in the CMV case. Our results confirm Conjec- 
tures 12.2.3 and 12.2.4 of [103J as well as Conjectures 12.2.5 and 12.2.6 in the 
(generic) special case that all gaps are open. 

For the case of OPUC with a single gap, the analog of Theorem 11.21 is known 
and motivated Simon's conjectures in [103| . In that case, the isospectral tori are 
labelled by a g (0, 1) and consist of {a^'^^ : A e 9ID)} where al^-* — Xa. Then 
dmi<^,T) ^ is equivalent to |a„| a and a„+i/a„ 1. This is often called the 
Lopez condition. Bello-Lopez [TU] proved the OPUC analog of Theorem 1 1.2 1 for this 
case if CTcss(>/) = cr(Jo) is strengthened to cr(J) — cr(Jo) (the analog of Rakhmanov's 
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result). The full analog for this special case appears in Simon (103j, Alfaro et al. 
[2], and Barrios et al. [9]. 

Associated to each two-sided p-periodic Jacobi matrix, Jq, is a polynomial, Aj^, 
of degree p, known as the discriminant. This is a classical object described in detail 
in the next section. It is usually defined as the trace of the one-period transfer 
matrix. It is also the unique polynomial (with positive leading coefficient) such 
that 

a{Jo) = {x : Aj„{x) e [-2,2]} 

In particular, two sequences of coefficients are isospectral if and only if they give 
rise to the same discriminant. 

The key to the proofs of our results is what we call the magic formula. Let J be 
a two-sided Jacobi matrix, then 

Aj„{J) ^ SP + S-P (1.35) 

if and only if J S 7}^. Here S is the right shift (cf. (|2.31|) ). In particular, ()1.35|) 
already implies that J is periodic! In the OPUC case, A is a polynomial in z and 
z~^. It turns out that A(C) is always self-adjoint; moreover, 

Ac,{C) = SP + S-P (1.36) 

if and only if C € Tcg . 

It has been previously noted that for periodic Jo, one has 

AM) = sp + s-p 

That this holds for some polynomial in Jq is in Naiman [751 [HH] • That the poly- 
nomial is the discriminant was found by Sebbar-Falliero [97]. After learning of 
our results, L. Golinskii has kindly pointed out to us that Naiman f80^ also has 
a theorem which implies any J obeying (|1.35|1 is periodic, the core of proving the 
converse. We will discuss this further in Section [S] 

Nonetheless, the two facts that make (|1.35p magical to us — that it characterizes 
the isospectral torus and that it is ideal for the study of perturbations — seem to 
have escaped prior notice. 

While J may be tridiagonal and C five-diagonal, both A(J) and A(C) are 2p+ 1- 
diagonal, that is, vanishing except for the main diagonal and p diagonals above and 
below. Thus, both will be tridiagonal if written as p x p blocks. The key to our 
proofs will be to extend results from the a„ = 1, 6„ = case to block tridiagonal 
matrices, and then use (|1.35p or (|1.36p to study perturbations of the periodic case. 

The magic formula is very powerful and opens up many other avenues for study: 

(a) Szego and Jost asymptotics for periodic perturbations and, in particular, the 
analogs of Damanik-Simon ^22^ . 

(b) Periodic analogs of the results of Nevai-Totik [SI] and its various recent ex- 
tensions f23l [TonfT05] . 

(c) Analogs of the Strong Szego Theorem for the periodic case following Ryckman's 
paper |95] for the Jacobi case. 

We should point out a major limitation of our results. If i? is a disjoint finite 
union of closed intervals in M (or 9D), one can construct an isospectral torus of 
Jacobi (or CMV) matrices whose recursion coefficients are almost periodic. As 
discussed in Section|2l these are strictly periodic if and only if the harmonic measure 
of each interval is rational. There are obvious potential extensions of Theorems ll.il - 
11.51 to this setting, but except for Theorem 11.11 (where the method of [69] applies) 
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and Theorem 11.21 (where Section [9] has some extensions) , we do not know how to 
prove them (or even if they are true) . There is no analog of A in the almost periodic 
case, so our method does not work directly. 

Here is the plan of this paper. Section [2] reviews the theory of the (unperturbed) 
periodic problem. In Section [3l we prove the magic formula for OPRL, and in 
Section m the magic formula for OPUC. While we will not discuss Schrodinger op- 
erators in detail here, we discuss the magic formula for such operators in Section [5l 
As we have mentioned, the magic formula brings block Jacobi matrices into play, so 
Section [S] discusses matrix- valued OPRL and OPUC — mainly setting up notation. 
Section[7|uses known results on Rakhmanov's theorem for matrix-valued orthogonal 
polynomials to prove a Denisov-type extension which we use in Section [8] to prove 
Theorem II. 2| the section also proves half of Theorem 11.11 Section [9] provides two 
results that go beyond the periodic case to prove Denisov-Rakhmanov-type theo- 
rems for special almost periodic situations. Section [TOl following |100j . proves the 
P2 sum rule of Killip-Simon 6IJ and the Ci sum rule for matrix-valued measures, 
and Section [11] uses these results to prove Theorems 1 1 . 31 and 11.41 Section[T2l explores 
what we can say if gaps are closed. Section [T3l proves analogs to the Lieb-Thirring 
bounds of Hundertmark-Simon '55] as preparation for proving Theorem ll.Sl in Sec- 
tion [TD Finally, Section [T5] discusses the OPUC results. 

Acknowledgements. It is a pleasure to thank Leonid Golinskii, Irina Nenciu, 
Leonid Pastur, and Peter Yuditskii for useful discussions. 

Note Added August, 2008. During the refereeing of this paper, Remling (in 
[93j ) ■ motivated in part by this paper, found a positive resolution of the conjecture 
that, in the language of our Theorem 19. 5( every set in is a Denisov-Rahkmanov 
set. His analysis depends on a very interesting theorem on right limits of Jacobi 
matrices with absolutely continuous spectrum - it provides a new approach to 
Denisov-Rahkmanov theorems. 

2. Review of the Periodic Problem 

In this section, we collect some of the major elements in the strictly periodic case. 
As this is textbook material, we forgo proofs and historical discussion. Full details 
can be found, for example, in [20l HQI [HI [1031 IMl [IIIl [IE] and the references 
therein. 

To discuss the strictly periodic case, we need to extend our operators to be two- 
sided, that is, to act on £'^{Z). In the Jacobi/OPRL case, we simply continue the 
tridiagonal pattern with parameters {a„,6„}„gz- Two-sided (or extended) CMV 
matrices are formed as C = CM, where C and M are doubly infinite direct sums 

£ = ...© e_2(a-2) © Oo(ao) ® 02(^2) e • • • (2.1) 
•••®e_i(a_i)®ei(ai)©--- (2.2) 

that are misaligned by one row/column, just as in (|1.5p - (|1.8p . 

We adopt the convention of indexing the elements of matrices so that 

Jii = bi Coo = ao Moo = -a-i (2.3) 

except A^oo = 1 in the one-sided case. 
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2.1. Transfer Matrices. Let J be a two-sided Jacobi matrix. A sequence {un} 
obeys (J — x)u = if and only if 

anUn+l + [hn - x)un + a„_iu„_i = (2.4) 

or, what is equivalent, 

'^-+^\^kJ (2.5) 



K{x) = -(''J- (2.6) 



with 

A \ = — P 

a„ \ a„ 

Note that the desire to have A„ depend only on one pair (a„, 6„) and have determi- 
nant equal to one resulted in the factors a„ and a„_i appearing in (|2.5p . (The same 
price is usually paid when writing Sturm-Liouville equations as first-order systems.) 
The choice is not the most common one (although it is used in Pastur-Figotin [86j). 
but we feel it is the 'right' one since, in particular, det(A„(a:)) — 1. 
In the OPUC case we define 

which encodes the recurrence relation (11.21): 



^"+^(^^^-Af„(z)f^:!f) (2.8) 



We will now explain the link to formal (i.e., not necessarily £^) eigenvectors of a 
two-sided CMV matrix, C. It is not as simple as (12.51) . 



Lemma 2.1. Suppose (C — z)u = with z ^ and let v = Z ^Mu where Z 
denotes the diagonal matrix with entries 



j odd 
j even 



and M is as in (12.21) . Then 



=zM2„+i(z)r^"+M =ilW(z)M2„(^)r'") (2.9) 

V2n+2 J \U2n+l I \V2r. ' 



Proof. The key observation used to verify (|2.9p is 

:;) = e(«„) ^ (;;)^M„(.)(^) (2.10) 

This follows by simple algebraic manipulations: 

CtnX + PnX' = Zy and PnX — OLnX' = y' 

x' = p^^{zy ~ Oinx) and y' = pnX — anx' 
x' ^ p~'^ {zy - Unx) and y' = p:^'^ {-a„zy + x) 
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With ()2.10|) now in hand, we may argue as fohows: 

(C ~ z)u^O 
'^F=^ zu = CMu 

V Z^^A4u obeys zu — CZv 
<F=^ Zv = Mu and zu = CZv 



f ZV2n-l 



= B(a2„-i 

= M2„-l(z) 
= M2„+l{z) 



U2n-1 
U2n 



U2n 
V2n 

U2n+1 
V2n+2 

which are the two parts of (|2.9 



V2n-1 



V2n+l 
U2n+l 



and 
and 
and 



ZU2n 
ZU2n+l 

f ZV2n+l 
\zU2n+l 

U2n+1 



= 0(a2«) 

= M2n{z) 
= M2„(z) 



V2n 
ZV2n+l 



U2n 
V2n 

U2r 

V2r, 



□ 



T(z) = 



(2.11) 
(2.12) 



2.2. The Discriminant. As in Sturm-Liouville theory, the discriminant is defined 
as the trace of the one-period transfer matrix: 

A(z) = Tr(T(z)) 

where 

' Ap(z)...A2(z)Ai(z) (OPRL) 
z-p/2Mp_i(z) • • • Mi{z)Mo{z) (OPUC) 

In the OPUC case, p is even. Also, the factor z^^/^ is there to cancel the extra 
factor of z on the left-hand side of (|2.9p . From a strictly OPUC point of view, 
it is more natural to omit this factor (as in |102[ I103j ): however, as the magic 
formula is an operator identity, we have elected to use the definition best adapted 
to this perspective. The only negative side effect of this choice is that our Lyapunov 
exponent (defined below) differs by — ^ log \z\ from that in jl021 1103] . 

For OPRL, the discriminant is a real polynomial of degree p with leading behavior 

p 



A{x) = (ai . . .Op) ^ 



l[ix-b,) + 0{xP-') 



(ai 



i=i 



rP-l 



0{X' 



For OPUC, it is a Laurent polynomial of total degree p with 



A(z) = A(l/z) 
so A is real- valued on 9D. Moreover, 

A(2) = (popi • • • Pp-ir'iz"^' + ■■■ + z-^l^) 



(2.13) 
(2.14) 

(2.15) 
(2.16) 



2.3. The Lyapunov Exponent. On an exponential scale, the behavior of formal 
eigenfunctions is determined by the Lyapunov exponent 



7(z)- lim — log||T"(z)|| 

ri^oo np 

= i log (spectral radius of T^z)) 
P 



(2.17) 
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= - log|A+(z)| 



(2.18) 



where X± are the eigenvalues of T{z) with the convention |A+| > |A_|. 
As det(r(2:)) = 1, these eigenvalues are the roots of 

- A(z)A + 1 = 



(2.19) 



which implies 



_ A(z)^ v/A^(z)-4 



(2.20) 




7(z) = ilog iA(z) + iVA(z)2-4 



(2.21) 



2.4. Gaps and Bands. Our recurrence relations admit bounded solutions for a 
given z if and only if A(z) e [—2, 2]. In the Jacobi/OPRL case, this is a collection of 
intervals in R. For CMV/OPUC, it is a collection of arcs in 9D. In either case, one 
may partition this set into p bands. These are the closures of the (disjoint) regions 
where A(z) € (—2,2). These can only intersect at the 'band edges', A~^({— 2,2}). 

The open gaps are the intervals/arcs that are complementary to the bands — 
excluding the two semi-infinite intervals in the OPRL case. When two bands touch, 
we refer to the common band edge as a closed gap. 

A^ — 4 has simple zeros at the edges of the open gaps and double zeros at the 
closed gaps; indeed, this is a complete list of its zeros. It is possible to distinguish 
whether these zeros correspond to A(z) — ±2 from the fact that there must be two 
zeros of A ± 2 between consecutive zeros of A =F 2 and the fact that A has positive 
leading coefficient. 

2.5. Spectrum. In both cases, the spectrum of the two-sided operator (acting 
on ^^(Z)) is the union of the bands: cr = A^^([— 2,2]). It is purely absolutely 
continuous and of multiplicity two. 

The spectrum of a two-sided p-periodic operator uniquely determines its dis- 
criminant; see Lemma 13.31 One consequence of this was noted already in the 
introduction: isospectral tori are the classes of p-periodic recurrence coefficients 
that lead to the same discriminant. 

In the case of a one-sided operator, the essential spectrum remains A~-'^([— 2, 2]); 
it is absolutely continuous with multiplicity one. In addition, up to one eigenvalue 
may appear in each open gap. 

2.6. Potential Theory. From the way they are defined, one can see that 7(2;) 
vanishes on the bands and is both positive and harmonic in the complement (in 
the OPUC case one must also exclude z = 0). This leads to the solution of the 
Dirichlet problem for a charge at infinity. 




(OPRL) 



(OPUC) 



(2.22) 



and so to the logarithmic capacity of the spectrum. 




(OPRL) 



(OPUC) 



(2.23) 
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2.7. Harmonic Measure. Taking normal derivatives in (j2.22|) leads to a formula 
for harmonic measure on a (aka equilibrium measure for the logarithmic potential), 



A'(x)| dx 



where supp((ii/) = cr = {z : |A(z)| < 2}. {Note: [103] has 1/p rather than 2/p, but 
that is an error.) 
Recognizing 

v/4- A2(a;) dx 

we see that the harmonic measure of each band is exactly 1/p. In particular, the 
connected components of the union of the bands all have rational harmonic mea- 
sure. This gives strong restrictions on sets that can be bands. In the OPRL case, 
rational harmonic measure of connected components is also sufhcient for a set to be 
the spectrum of a periodic Jacobi matrix. In the OPUC case, there is an additional 
condition needed: after breaking the bands into arcs of harmonic measure l/p, the 
harmonic midpoints of these intervals must obey 11^=1 Cj — 1- Clearly, the 

condition on the harmonic midpoints can be achieved by simply rotating a. Dis- 
carding this condition gives rise to Verblunsky coefficients that are p-automorphic, 
a„_i_p = e*'^a„, rather than p-periodic. 

2.8. Thouless Formula. Harmonic measure appears naturally in the theory in 
several other ways. It is the density of states measure: 



if (OPRL) 



^ ' - arccos — ^ (2.25) 



1 ^ r 

i™^2iV ^ /(•^)"" = / f(^)dHx) (2-26) 

for every polynomial (or continuous function) /. The same formula holds with C 
replacing J. This connection, or more precisely the resulting expression for 7(2;) 
in terms of Cap(CT) and the logarithmic potential of dv^ is known as the Thouless 
formula. 

Two further characterizations of dv involve the orthogonal polynomials, dv is 
the weak limit of j^^n=o P^^^) d^i.^) (resp. Y.nZo\Vn[e'')?dti[e)). It is also 
the limiting density of zeros, that is, the weak limit of the probability measures, 
dvm which give weight 1/n to each of the zeros of p„ (resp. ipn)- These two charac- 
terizations are closely linked to (|2.26p : however, in the OPUC case one should keep 
in mind that for each n, the zeros of ipn lie strictly inside D. 

2.9. Floquet Theory. Looking at the eigenvalues of T, one sees that when A+ 7^ 
A_, there is a basis of formal (i.e., non-^^) eigenfunctions obeying 

Um+kp — ^±Um (2.27) 

If A+ = A_, which happens precisely at the band edges, then both are ±1. If the 
edge abuts an open gap, there is only one eigenfunction obeying (|2.27p since T has 
a Jordan block structure. At closed gaps, T = ±1 and so all solutions obey (|2.27p . 
Solutions obeying 

Um+kp = e'^^Um (2.28) 
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are called Floquet solutions and e*^ is called the Floquet index; they have much 
the same role as plane waves in Fourier analysis. Since A_ = A^^, if (|2.28|1 has a 
solution, then e~'^ is also a Floquet index. 
In the OPRL case, 

([2:281) holds ^=> A{x) = 2 cos 6* (2.29) 

Thus, by the discussion above, for each 6 G (0, tt), (|2.28p or (|2.29p holds for exactly 
p values of x: xi{9) < X2{0) < ■ ■ ■ < Xp{9). These Xj{d) are known as the band 
functions. 

The changes in the OPUC case are purely notational: 

(|2:28l) holds A(z) = 2 cos 6* (2.30) 

For 9 G (0, tt), this has p solutions all of which lie in 9D. 

2.10. Direct Integrals. Let S denote the right shift, 

{Su)n = Un-1 (2.31) 

If the sequences of coefficients are p-periodic, then J (or C) commutes with S^, 
which means that the two operators can be 'simultaneously diagonalized'. We 
elaborate this point in the OPRL case; the OPUC is almost identical. 
Let us write 

nr := |%^£ = i^([0,27r),|i;C^) 
where ig is the p-dimensional Hilbert space 

p 

£l = {u\Un+p = e'^U„} {u\v)g = ^ UnV„ 

From Fourier analysis, there is a unitary operator T: £^(Z) Tip so that 
TS'^T^^ is multiplication by and TJT~^ acts fiber-wise (i.e., on each £g) 
as a p X p matrix, J (9). In particular, the eigenvalues of J (9) are the solutions of 
(|2.29p . that is, they are the band functions Xj{9). 

2.11. Hyperelliptic Riemann Surfaces. As \/A^ — 4 appears repeatedly in the 
theory, it is natural that the associated Riemann surface should enter the analysis. 

— 4 has simple zeros at the edges of open gaps and at inf cr(J) and sup(T(J). 
It has double zeros at the closed gaps. Let i denote the number of open gaps, 
then \/A^ — 4 has square root singularities at 2(^+1) points, and so its natural 
analyticity domain is the genus £ Riemann surface, S, obtained by taking two 
copies of C \ (t{J), gluing at the bands and adding points at oo. There is a natural 
projection tt: S CU{cxd} which is 2 to 1 except at the branch points of VA^ — 4. 
A similar analysis works for OPUC, but now there are £ gaps and the genus is ^ — 1. 

2.12. Minimal Herglotz and Caratheodory Functions. For a half-line peri- 
odic Jacobi matrix, the m function is defined by 

m{z) = {So,{J - z)-'^6o) VzeC\cr(J) (2.32) 

This can be shown to obey a quadratic equation with polynomial coefficients 

A{z)m{zf + B{z)m{z) + C{z) = (2.33) 



14 



D. DAMANIK, R. KILLIP, AND B. SIMON 



Moreover, these coefRcients can be chosen to obey 

_ \AC = - 4 (2.34) 

This imphes that m(z) has meromorphic continuation to S. Indeed, m has minimal 
degree (i.e., degree £ + 1 in the OPRL case and I in the OPUC case) among aU 
meromorphic functions on S that are not of the form g o tt with g meromorphic 
on the Riemann sphere. It can be shown that there is a one-one correspondence 
between minimal meromorphic functions obeying ImTO(z) > if Imz > and 
m{z) = —z~^ + 0{z~'^) on the top sheet of S and all periodic Jacobi parameters 
with the same A. (We call these minimal Herglotz functions.) 
There is a similar description for OPUC, but now one uses 

F{z)^{SoAC + z){C-z)-^So) (2.35) 

which obeys Re_F(2) > if |z| < 1 and F{Q) — 1. Again F obeys a quadratic 
equation, showing that F has a meromorphic continuation to S of minimal degree, 
and again there is a one-one correspondence between all {q;„}^^q with the same A 
and all minimal Caratheodory functions. 

2.13. Dirichlet Data. One can describe the set of minimal Herglotz functions in 
terms of their poles. For each open gap, {Gj}j^i, ^^^{Gj) = Tj is a circle since tt 
is 2 to 1 on Gj and one-one on Gj \ Gj. A meromorphic Herglotz function has £+1 
simple poles, one at oo on the second sheet and the other £, one in each Tj. Thus, 
the set of meromorphic Herglotz functions is homeomorphic to X^^^^ Tj under the 
bijective map from such functions to its poles. A similar analysis holds for OPUC 
but now there is no pole at infinity, there are £ gaps, and Xj^^Tj describes the 
possible poles. The difference is that for OPRL, the dimension of the torus is £, 
and for OPUC it is ^- 1. 

2.14. Isospectral Tori. By combining the bijective maps from periodic OPRL to 
minimal Herglotz functions and of such functions to Dirichlet data, we see for a A 
of period p with £ gaps, 

{{an,bn)n^i ■ the discriminant is A} 

is an ^-dimensional torus in M.'^p. Generically, £ = p — 1. In the OPUC case, 
generically £ = p and the torus is naturally embedded in C^. This torus is the 
isospectral torus which we will denote by T or Tj^ if a given periodic Jq underlies 
our construction. For clarity of exposition, we will typically blur the distinction 
between p-tuples (a„,6„)^^]^ and the corresponding infinite sequences {a„,5„}„gz 
of period p. Because of our perturbation theory viewpoint, we use Jo to label the 
torus, but we emphasize that from another point of view, the torus is associated to 
the set iToss(Jo) and not to Jq. 

2.15. Isospectral Flows. The fact that spaces of p-periodic coefficients foliate 
into tori suggests that there is some kind of completely integrable system in the 
background. That is true: it is the Toda flow in the OPRL case and the defocusing 
Ablowitz--Ladik flow in the OPUC case. Since we will not need these below, we say 
no more about them, but see Chapter 6 of [106' for the OPRL case and Section 11.11 
of [103] for OPUC. 
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3. The Magic Formula for Jacobi Matrices 

Our goal in this section is to prove 

Theorem 3.1. Let Jq be a two-sided p-periodic Jacobi matrix with discriminant 
Ajg (x) and isospectral torus Tj^ . Let J be a two-sided [not a priori periodic) Jacobi 
matrix. Then 

Aj„(J) ^SP + S-P ^ J e Tj, (3.1) 
where S is the right shift, (j2.3ip . on ^^(Z). 

We provide two proofs of the 'harder' direction or rather of 

Aj„(J) = 5*^ + 5-^^ J is periodic (3.2) 

which is the key step. Our first proof is immediately below; the second, suggested 
to us by L. Golinskii, appears after Lemma 13.41 

Lemma 3.2. Let t ^ 1,2, ... . Then 

{0 k>t 
O'mO'm+l ■ ■ ■ 0,m+k-l k=£ (3-3) 

amOm+l • ■ ■ Clm+e-2{bm + bm+1 + ' ' ' + bm+l-l) k — £ — 1 



Proof. Writing 



(^J )m,m4-fc — Jm,,iiJii,i2 ' ' ' Jiit-i.m+k (3-4) 



We see that since J is tridiagonal, all terms are zero if fc > ^, that we must have 
(with io = m, ig — m -\- k) that iq — iq-i = 1 for q = 1, . . . ,£ \i k = £, and that if 
k ^ £ — \, iq — iq-i = 1 for all but one q € {1, ...,£} and it is zero for that q. □ 

Lemma 3.3. If J and Jq are periodic, then a{J) = (j{Jo) if and only if Aj — Aj^ . 

Remark. This lemma says that the spectrum determines the discriminant and vice 
versa. That the discriminant determines the spectrum is elementary: a = {x : 
A{x) £ [—2, 2]}. Therefore we only prove the other direction — indeed, we give two 
proofs. 

First Proof. Harmonic measure dv is intrinsic to the set a; it is the solution of an 
electrostatic problem there. But then dv determines A via ()2.24|) . □ 

Second Proof, a determines the gaps — even closed gaps — via harmonic measure. 
The gap edges determine the zeros of A — 2 and so A — 2 up to a constant. The 
zeros of A + 2 then determine the constant. □ 

Proof of Theorem V3.1\ For all Q £ [0, 2tt), J (9) is self-adjoint and so diagonalizable. 
Moreover, the eigenvalues of Jo{6) are precisely the roots of A(a;) = 2cos(6'). Thus 

Aj„(Jo(^)) = (2cos0)l (3.5) 

But then Aj„(Jo) and S'' + 3^^ both have direct integral decomposition with fibers 
(2 cos 6*) 1, so 

AM)^SP + S'P (3.6) 
Since J £ Tjg ^ Aj = Aj^, this proves <= in (|3.ip . 

Now suppose LHS of holds. By ((2T4| . and 

{SP + S-PU,ra+p - 1 iS" + S-P)m,m+p-l = (3.7) 
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this implies 
and 

In particular, 
which lead to 

dm — Qm+p — ^m+p (3.10) 

SO J is periodic. 

Since J is periodic, Aj(J) ^ S'p + S~P; moreover, Aj„(J) = 5^ + S'P by 
hypothesis. Thus we learn that applying the polynomial Aj — Aj„ to J gives 
zero. By the k — £ case of p.3p . it must therefore be the zero polynomial, that is, 
A J = Ajg. Lemma 13.31 now completes the proof. □ 



•••a™+p-i (3.8) 

p—l p—1 

E^™+^ = Ed (3-9) 

3=0 j=o 



p-i 

3=0 



Remarks. 1. Showing that J was periodic only required equality in Ajq(T) — 
SP + S^P, for the two most extreme upper (or lower) diagonals. Nevertheless, 
J G X/o requires equality everywhere. 

2. We need not suppose a priori that each On > and can allow some On = 
(J can still be defined on ^^(Z)), for dSj]) implies that if LHS of ([33]) holds, then 
each a„ > 0. 



We now turn to our second proof of (|3.2p . 

Lemma 3.4 (Naiman [80]). Lei A be a two-sided (bounded) infinite matrix of finite 
width {i.e., for some w, we have that \k — £\ > w ^ A^i = 0). Suppose 

[A,SP + S-P] = Q (3.11) 

for some p, then 

[A,SP]=Q (3.12) 

Remarks. 1. This is Lemma 2 in [80 ; no proof is given. 

2. [A, B] = AB-BA 

3. (|3.12p has an equivalent form: 

[A, SP] = Q ^ Ak+p.i+p = Ak,i for aU fc, t (3.13) 

4. As (|3.13p shows, [J, S^] — for a Jacobi matrix if and only if and hk are 
p-periodic. 

Proof. Since A has finite width, we can find diagonal matrices , -Dfej+i, . . . , Dk^ 
with Dfej 9^ ^ ZJfej , so that 

A^Y^ DjS' (3.14) 
Since Dj is diagonal, so is S^DjS^P. Thus 

{SP + S-P)A = J2 {SPDjS-P)S^+P + J2 {S-''DjSP)S^-P 
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^2 ^2 

j=ki j = ki 

Since the composition (I3.14|l uniquely determines each Dj, implies 

SPDk,S-P = Dk, 

that is, 15^2 is periodic. Thus, D^^Sp commutes with Sp + S~p^ so we can remove it 
from (|3.14p without losing p. lip . This shows inductively that each Dj is periodic. 

□ 

Second proof of (13. 2p . J commutes with A(J), so (|3.12p holds. □ 

Our next goal is to compare dm{{a, b), T) given by (|1.34p and dm{{a, b), T) given 
by (|1.18p . As well as satisfying natural curiosity, this relation also plays an impor- 



tant role (via Theorem I11.13P in the proofs of Theorems 11.31 and 11.41 

To capture the essence of what follows, let us pause to ponder the following: 
suppose dm{{a, 6), T) = for all to, does this mean that (a, b) £ T? The hypothesis 
tells us that each length-p block belongs to the isospectral torus; it does not a 
priori even guarantee that the coefficients are periodic. Example 14.51 shows that 
periodicity can fail in the OPUC case. However, such problems do not arise for 
OPRL. The reason is simple: within the isospectral torus, ai, . . . ,ap_i determines 
Op and 6i, . . . , 6p_i determines bp. 

Proposition 3.5. Given a p-periodic Jacobi matrix Jo, ^ < q < oo, and e > 0, 

there is a constant C so that 

ei-P||J„,((a,6),TjJ||^,< ((a, 6), T/J ||^,< C||d„, ((a, 6), TjJ ||^, (3.15) 

for all sequences {(a„,6„)} obeying > a„ > e > 0. All I'' norms are over 
me {1,2,3,...}. 

The key input is 

Lemma 3.6. Given {(a„, 6„)} obeying > a„ > e > 0, 

n—p+l 
r—m 

for all n > m. The constant C depends only on e. 

Proof. The proof is by induction on n. For TO<n<TO+p — 1, the result is 
immediate from the definition of d,n . 

For n > m + p — 1, we consider the functions 

/(ai, ...,ap):=j2 H{a,) - log(af )] g{b,, . . . , 6,) := [b, - bf^] 



These vanish on Tjg , as explained in the proof of Theorem 13.11 
As g is Lipschitz (with constant 1), 

\bn — bn-p\ = \g{bm ■ . . , bn-p+l) — .g(&„-l, . . . , bn-p)\ 

< \g{bn,...,bn-p+i)\ + |.9(6„_i, . . . ,5„_p)| 

< d„_p+i((a,6),T/o) + dn-p{{a,b),Tjg) 
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In a similar way, 

|log[a„] - log[a„_p]| = |/(a„, . . . ,a„_p+i) - /(a„_i, . . . , a„_p)| 

leads to 

((a,5),T/(,) + d„_p((a,6),TjJ] 
Combining these two inequalities gives 
L _a(o)| I |5 _;,(0)|<|a _ | , _ ^(o) i 

+ (1 + Ce)[d„_p+i((a, 6),T7„) + dn-p{{a,b),Tjg)] 
which completes the proof of the inductive step. □ 

Proof of Proposition \3. 51 The left-hand inequality in (|3.15p follows immediately 
from the definitions of dm and dm] we focus on the second inequality. 

Choose (a(o),5(°)) minimizing dm{{o,,b),Tjg); strictly, this amounts to a (incon- 
sequential) change in Jq. Applying Lemma 13.61 in the definition of d,„ gives 

oc m-\-k 

dm{{a,b),T.j„) < ^dm{{a,b),Tjo)+CY, E dr{{a,b),Tj„) 

k—0 r—m 

oo 

<C'J2e-^dm+,{{a,b),Tj,) 

The proposition follows because convolution with e'^-' X[o.oc){j) is a bounded oper- 
ator on all spaces. □ 

4. The Magic Formula for CMV Matrices 

Our goal in this section is to prove 

Theorem 4.1. Let p be even and let Co be a two-sided p-periodic CMV matrix with 
discriminant Acoiz) and isospectral torus Tca- Given a two-sided {not a priori 
periodic) CMV matrix, C , 

Aco (C) ^S^ + S-P ^ Ce Tc, (4.1) 

Remarks. 1. Notice that since C is unitary and A(e*^) is real, Aco(C) is self-adjoint. 

2. By (|2.16p and the fact that C is five-diagonal, Ac(,(C) has 2(p/2) diagonals 
above/below the main diagonal. 

3. As in Section [3l we will first present our initial proof that 

Aco (C) ^SP + S^P => {an} is periodic (4.2) 
and then a proof based on Golinskii's suggestion. 
Lemma 4.2. We have: 

{C')m.m+k - iC-')m,m+k =0 if k > 21 (4.3) 

(CO 2m, 2m +2£ — P2mP27n + l ■ ■ ■ P2m+2l-l (4.4) 

(C^)2m+l,2m+2£+l = (4.5) 

{C'%m,2n^+2^ = (4.6) 

{C ^)2m+l,2m+£+l = P2m+lP2m+2 • • • P2m+2l (4-7) 

(CO 2m,2m-|-2£-l — P2mP2m+l ■ ■ ■ P2m+2l-2a2m+2l-l (4.8) 
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(C^)2m+l,2m+2£ — —a2mP2m+l ■ ■ ■ P2m+2l-l (4-9) 
2 m, 2m +2£-l — —Cl2m,-lP2mP2m+l ■ ■ ■ P2m+2£-2 (4-10) 
2m+l,2m +2i = P2m+1 • • ■ P2m+2e-ia2m+2e (4-11) 

Proof. As C and Ai are tridagonal, is a product of 2£ tridiagonal matrices, so 
(|4.3p is immediate. 

We will prove the results for C^. The results for are similar if we note 

e(a)-i = e(a) (4.12) 

since Q is unitary and symmetric. 
Equation (|4.4|) follows from 

^2m,2m+l = P2m ■A^2m+l,2m+2 — P2m+1 (4-13) 

and (|4.5p from 

>C2m+l,2m+2 = (4.14) 

Because of (|4.14l) , the only way for to get from 2m to 2m + 2^ — 1 is to increase 
index in the first 2^—1 factors, which leads to (|4.8p . For the same reason, to get 
from 2m + 1 to 2m + 2£, the last 2£ — 1 factor must increase index, leading to 
gH). □ 

Lemma 4.3. IfC and Co are p- periodic, then<j{C) = cr(Co) if and only if l^c — ^Co- 

Proof. Either proof of Lemma 13.31 carries over with no change. □ 

Proof of Theorem The proof that 

^cACo) ^ + S-P (4.15) 

is identical to the proof of ()3.5p . 
For the converse, suppose 

AcoiC) ^ SP + S-P (4.16) 

In particular, 

Aco(C)2m,2m+p-l = (4.17) 

By (|2.16p and Lemma 14.21 this implies (recall p is even) 

2 m Pm+l ■ ■ ■ P2m+p-2){a2m+p-l - a2m~l) = 

SO 

a2m+p-l = a2m~l (4-18) 

Similarly, since 

^Co{C)2m+l,2m+p — 

we get 

2m-\-l • ■ • P2m+p-l){a2m+2 ~ a2m+p-l) = 

which leads to 

a2m+p = Q-im (4-19) 

Thus, a has period p. That C G follows from Lemma 14.31 and the same 
argument used in the OPRL case. □ 

Next, we give a proof using Naiman's lemma. We will need 
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Lemma 4.4. Let C be the extended CM V matrix associated to {anjJ^L-oo- Let p 
be even. If S^C — CS^, then 

an+p = an (4.20) 

for all n. 

Proof We have that C^. ^j+i + C^j 2j+2 = Plj (see (4.2.14) of [TU^ ). so p2j is 
periodic. Thus, C2j 2j+2/p2j ~ P2j+i is also periodic. So a2j+i = C2j 2j+i/p2j is 
periodic as is a2j = C2j+i 2j+2/(-p2j+i)- □ 



Second proof that (|4.2p holds. C commutes with + S~p, so by Naiman's lemma 
(Lemma 13. 4p . which did not require that A be self-adjoint, S^C = CS^, which 
implies a is periodic by Lemma 14.41 □ 



We now turn to the OPUC version of Proposition 13. 51 As noted in the introduc- 
tion, it is not sufficient to sum over exactly one period: 

Example 4.5. (0,^,0,^,0,...) and (0, —^,0, — ^, 0, . . . ) are in the same isospec- 
tral torus, namely, the one with p ~ 2 and 



A(z) = ^|(z + z-0 

Now consider a — (0, i, 0, — i, 0, i, 0, — i, . . . ). If d„i{a, Tc„) were defined as sum 
from fc = to p — 1, it would be zero for all m, but dm{oL, Tcg) is not small. □ 

The problem, as this example shows, is that for sequences in Tc,,, {ckq, . . . , ap^2) 
does not determine ccp-i. But by periodicity, Uq, . . . , ctp-i determines Up. Thus, if 
we define 

p 

J„i(a, a') := ^ \a„r+k - a'^+k\ (4-21) 
k=0 

then 

|am+p - ami < dm{a,Tca) 
Plugging this into the proofs of Lemma [3.61 and Proposition 13.51 leads quickly to 

Proposition 4.6. Let Co be a fixed periodic CMV matrix, then 

e~P\\djn{a,Tc„)\\^^< \\djn{a,Tc„)\\f^< C\\dra{a,Tc„)\\f,^ (4.22) 

for any sequence of Verblunsky coefficients {a„}. 

5. The Magic Formula for Schrodinger Operators 

In this section, we want to illustrate the potential applicability of our central 
idea to the spectral theory of one-dimensional Schrodinger operators, 

H = -— + V{x) (5.1) 
ax 

However, we will not pursue the applications in this paper. 

We wiU suppose V € iinifi that is, sup^ J^^^ \ V{y) \ dy < oo. In that case, V is 
a form bounded perturbation of on L^(K, dx) with relative bound zero, so H 

is a self-adjoint operator. Its form domain is the Sobolev space H^{M.). 

We need to say something about periodic Schrodinger operators. Suppose Vq 
has period L, that is, 

Vo{x + L) ^ Vo{x) (5.2) 
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For arbitrary V in Lj^^, and E E C, let ud{x, E;V) and un{x, E;V) (we will 
often drop the V if it is clear which V is intended) be the solutions of 

- u" + Vu^Eu (5.3) 

obeying the boundary conditions 

ud{0) = Q ua{0) = l un{0)^1 w^(0) = (5.4) 

There is a unique solution of (|5.3p in distributional sense which is absolutely con- 
tinuous. 

The transfer matrix that updates solutions of (|5.3p (with data written as (",)) 

is 

T(x,i?;l^)=fTh^^! r^^^O (5.5) 
\Un[x,E) u'jj{x,E)J 

det(T) = 1 by constancy of the Wronskian. For periodic Vb, we define the discrim- 
inant by 

Ava{E)=Tr{TiL,E;Vo)) 

= UN{L,E) + u'a{L,E) (5.6) 

As in the OPRL and OPUC cases, it is easy to see for the whole-line operator 
that 

a(-i^ + Fo)=A^„H[-2,2]) (5.7) 
and is purely absolutely continuous. Moreover (see, e.g., [92]), if 

{Syu){x) = u{x - y) (5.8) 

,2 

then H — —-^ + Vq commutes with Sl and so has a direct integral decomposition, 

r® d0 

H = j H{B)- (5.9) 

whose fibers, H{d), are the operator (|5.ip on [0, L] with 

u{L) = e'\{0) u'{L) = e^%'(0) (5.10) 

boundary conditions. El(Q) has purely discrete spectrum (i.e., {H{9) + i)^^ is 
compact); the eigenvalues are precisely the solutions of 

A(i;) = 2cos(6') (5.11) 

Two periodic potentials of period L are called isospectral if and only if they have 
the same A. As in the Jacobi and CMV cases, the spectrum determines A, but this 
is more difficult to prove in the Schrodinger case. It is also known ([78] for nice Vb's; 
[T51 BTl [551 [5T]) that the set of V's isospectral to Vq is a torus of dimension equal 
to the number of gaps which is typically infinite, so we will refer to an isospectral 
torus, Tv„- We can now state the main result in this section: 

Theorem 5.1. Let Vq be periodic obeying (j5.2[) and let Ay^ be its discriminant 

1 

and its isospectral torus. Let V be in L^^^ ^^^^ on R and H = -j^ + V. Then 
Av„iH) = Sl + S-l ^ Ve (5.12) 
Here S±l denotes the shift operator, as in (|5.8p . 
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Remarks. 1. Avo(iJ) is defined by the functional calculus. 

2. As in the last two sections, we will provide our initial proof that 

Avo (H) ^ Sl + S-L^V periodic (5.13) 

and then a simpler proof using an analog of Naiman's lemma. This argument 
does not require Theorems 15.21 and 15 . 31 and the considerable machinery their proofs 
entail. That said, to show A Vo {H) — Sl + S-l plus V periodic implies V £ Ty^ 
does require Theorem 15. 3| but it should be noted that one can prove Theorem 15.31 
fairly easily without needing transformation formulae of Delsarte, Levitan, Gel'fand, 
Marchenko type. 

We need two preliminaries whose proofs we defer to later in the section. We first 
make a definition: 

Definition. For any y > 0, TZy consists of operators on L^(M) of the form 

{Af){x) ^lf{x + y) + y{x-y)+ f K{x, z)f{z) dz (5.14) 

where K is continuous and uniformly bounded on {{x, z) : \x — z\ < y}. 

Note. It can happen that K{x,x ± y) 7^ 0, so if we think of K as an integral 
kernel on K x R, it can be discontinuous at |a; — z| = y. 

Theorem 5.2. // Vq is L-periodic and V in L\^^ then ^A{H) e TZl and 

Kix, x + L) = -\ f ^ {V{z) - Vo{z)) dz (5.15) 

J X 

Note that (|5.15p describes the 'matrix elements' of Avo{H) — {Sl + S-l) that 
are farthest from the diagonal. Indeed, just as in the other cases, one does not 
need the full statement Ay,, [H] — Sl ~\- S-l to see that V is periodic, only that 
(/, (A(iJ) — Sl — S-L)g) = for / supported near xq and g near xq ^ L (for all 

Xq). 

Theorem 5.3. A{E) is an entire Junction which obeys 

(i) |A(£:)| < Cexp(Lv^) (5.16) 

(ii) lim =1 (5.17) 

Proof of Theorem \5.1\ If V G Tv;,) then Ay — Ay^, so for the -4= direction we need 
only prove 

Av,{£i + Vq)^Sl + S-l (5.18) 
As before, this is equivalent to Ay(,(iJ(0)) — 2 cos 6* which follows from (|5.1ip . 

Conversely, if Avo{H) = Sl + S-l, then from Theorem 15.21 and the periodicity 
of Vb, we see 

fX+L 

V{z) dz — constant (5.19) 

This implies that V{x + L) ~ V{x) = for a.e. x, that is, V is periodic. 

If H{9) are the fibers of H in the direct integral decomposition, Avo{H) = 
Sl + S^L implies 

Ayo(iJ(6i)) = 2cos6' (5.20) 
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SO, if A is the discriminant for V, we have A(z) — ±2 Avo{z) — ±2. Moreover, 
implies cr{H) C + y^)^ 

so any double zero of A ± 2 is a double zero 

of Ay;, ± 2. It follows that 

Al (z) - 4 

is analytic. 

Since Ay;, and A are entire functions of order ^ (by Theorem l5.3[) . g{z) is of the 
form 

g(z) = Cn(l-^) (5.22) 

where zi < Z2 < • • • are bounded from below. By (|5.17|) . ImiE^-oc giE) = 1, 
which implies g = 1, that is, A = Ay^. □ 

The argument used at the end of the proof to conclude that missing zeros cannot 
occur is reminiscent of ideas connected with the Hochstadt-Lieberman [S^ and 
related theorems [47l SS] . 

We now turn to the proofs of Theorems 15.21 and 15.31 A critical role is played by 
the wave equation and the transformation operator formalism of Gel'fand-Levitan, 
further important work is due to Delsarte, Levin, and Marchenko; see the book of 
Marchenko [76] for references and history. 

Define for s > 0, 

Cs{z) = cos(sVi) S,{z) = z-^/'^sm{sy/?i (5.23) 

which are entire functions of z bounded on (a, oo) for any a e M. Thus Cs{H) and 
Ss{H) are bounded operators for any H that is bounded from below. We will need 
to study the form of Cs{—-^ + V). For bounded continuous V, this is discussed 
in Marchenko 76j. While his proofs extend to the L^^^ case, it seems simpler to 
sketch the ideas: 

Proposition 5.4. C" := Cs{--^) € Us', indeed, 

{C^J){x)^\[f{x + s) + f{x-s)] (5.24) 

IfS^,:=Ss{~£,), then 

{S'lf){x) = \ f{y)dy (5.25) 

J X — s 

Remark. If w{x, s) :— {C^f){x) + {Ssg){x), then w obeys the wave equation (^ — 

■§^)w — with initial data w(x, 0) = / and dsw{x, 0) = g{x). Thus the proposition 
basically encodes d'Alembert's solution of the wave equation. From this point of 
view, Theorem 15.21 is connected to finite propagation speed for the wave equation. 

Proof. Since cos is even, 

cos(s|fc|) = cos(sfc) = i (e**-'" + e-'''') (5.26) 
()5.24p is just the Fourier transform of this. ()5.25|) follows from 



S,{z) = / Ct{z)dt (5.27) 

"'0 

and ([5:241) . □ 
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We are heading towards 

Theorem 5.5. Let V G Ll^uniii^) '^'^^ let H = --£:i + V. Then Cs{H) e TZs and 
the associated kernel of (|5.14[) obeys 

V{u)du (5.28) 



K,{x,x + s) = -\ j 

and for each t S (0, oo), 



sup \K,{x,y)\<(x, (5.29) 

x,y,\s\<.t 

In addition, 

{Ss{H)f){x)= f ^ Ls{x,y)f{y)dy (5.30) 
<y x—s 

where 

Ls{x,x + s)^\ (5.31) 

Lemma 5.6. It suffices to prove Theorem 1 5. 51 for s small. 

Proof. Since cos(2m) — 2cos^(u) — 1, one sees 

C2s{A) = 2Cs{Af - 1 (5.32) 

Thus, if Cs ^TZs, one sees C2s G 7^28 and 

K2s{x, y) = Ksix, y + s) + K^ix, y - s) + Ksix + s,y) + Ks{x - s, y) (5.33) 

where K{x, y) = Q ii \x — y\ > s. Thus 

K2s{x,x + 2s) = Ks{x,y + s) + Ks{x + s,y + 2s) (5.34) 

This shows that if the formula is known for \s\ < T, one gets it successively for 
2T, 4r, ST, ... . 

Using ([5^ . one sees that the result for Cs{H) implies ((OOl) and ((53T|) . □ 



Proof of Theorem 15.51 If A is a bounded self-adjoint operator on Ti which is bounded 
from below, and B is the operator on 7i 7i given by 



then 

" -^5* (A) Ct(A) 
This formula can be checked by showing that the right side of (|5.36p is a bounded 
semigroup whose derivative at i = is _B. DuHamel's formula for A, A bounded 
says that 

„tB _ tB I / „sB 



e — e 



f e'^{B-B)e^'-'^^ds (5.37) 
Jo 

= e*^+ / e"^(B-S)e(*-")^(is (5.38) 
Jo 

Using (|5.36p . we obtain 

CtiA) = Ct{A) - f Ss{A){A - A)Ct-siA) ds (5.39) 
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= Ct{A)^ [ S,{A){A- A)Ct^M)ds (5.40) 

By taking limits, it is easy to obtain these formulae for A — — A ~ + ^ 
with V bounded. By obtaining a priori bounds below depending only on certain 
norms of V, we get estimates for V in and so, using the lemma, prove the 
theorem. 

By iterating (j5.40p , one gets an expansion (which converges if V is bounded and 
whose estimates then extend), 

oo 

Ct{-£ + V{x))^ci"^+J2ci-^ (5.41) 

n=l 

= (-1)" / Si''^^VSi"J...VSiycE.^..._^y^i,---dsn (5.42) 

Jo<si-\ hs„<t 

Apply the integrand in C^"'' to a function / and evaluate at x for fixed si, . . . , s„. 
Each si^^'V evaluates V at points and integrals using ()5.25p . The integrands in V 
are in the interval {x — t,x + t), so if we take absolute values, we see this integrand 
is bounded by 

/ rx+t \ n 

U J ^ \Viy)\dy\ [lf{x + t~si sn) + yix-t + si + --- + s„)] 

Now we can do the integral over si, . . . , s„. For t — Si — • • • — s„ fixed, the new 
integrand is independent of Si, . . . , s„_i and is bounded by <"^^. We find 

\{d,"^f){x)\<r-'(U \V{y)\dy) / \f{y)\dy (5.43) 

\ Jx-t J Jx-t 

Moreover, C^"-* has a continuous integral kernel K["'\x,y) supported in \x — 
y\ < t. Since V is uniformly locally L^, by taking t small, we can be sure 
supj. \ {y) \ dy < I, which yields uniform convergence of K^"'' to a uniformly 

bounded kernel. 

By (|5.27p . we get (I5.30p from Cs(H) E TZs, and (|5.3ip comes from noting that 
\Ls{x,y)^^\<{.8-\x-y\) sup \K^ix,y)\ (5.44) 

x.y,u<s 

Finally, using (|5.39p . we see that 

Kt{x,x + t)——^ / Ls{x,x + s)V{x + s) ds 
Jo 

proving □ 

To complete the proofs of Theorems 15.21 and 15.31 (and so Theorem 15. ip , we need 
the transformation formulae of Delsarte, Levitan, Gel'fand, and Marchenko [7B]: 

Theorem 5.7. If V E L^([0, i?]) for R < oo, then there exist functions K^jK^, 
in {{y, x) : < y < X < R} so that for < x < R, 



un{x,E)^Cx{E)+ KN{y,x)Cy{E)dy (5.45) 
Unix, E) = Sx{E) + / Koiy, x)Sy{E) dy (5.46) 
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Moreover, 

Kd{x,x)^Kn{x,x)^\ ( V{t)dt (5.47) 

Remarks. 1. These formulae are in Marchenko [TBI P- 9 and (1.2.28)]. He supposes 
V is continuous, but his proof works if V is indeed, see Remark 2. 

2. Defining ux{x,k) — uxixjk"^) for X = D,N and Qx{x,y) as the Fourier 
transform of ux in fc, we see (I5.3P becomes 



with initial conditions 

Qn{x = 0, y) = Siy) Q'^ix ^0,y)^0 
Qd{x = 0, j/) = Q'd{x = 0, j/) = (5(2/) 

Thus, Theorem 15.71 is essentially Theorem 15 . 5 1 with a time-dependent V used. 
3. By (|5.6p . (|5.45p . and (|5.46|) . we obtain a critical representation for A: 

r-L 

\dt 



(5.49) 



A(S) = 2C,(i;) + / L,it)Ct{E)< 
Jo 

+ I L2{t)St{E)dt + Kd{L,L)Sl{E) 
Jo 

where Li,L2 are continuous in [0,-L]. Indeed, 

Li{t)=KN{t,L) L2{t) = ^Koit^x) 

Proof of Theorem \5.S[ The analyticity is immediate from (|5.49p as is (|5.16p given 

\C4E)\ + \S4E)\ < Cexp(a;^) 
Moreover, since for t < L, 

lini ^^=0 and lim = 

B^-oo Cl{E) E^-oo Cl{E) 

we have (|5.17p . □ 

Proof of Theorem\SIE By ([05)) and Theorem [531 | A(-^ + y) is in 7^t. More- 
over, the only terms contributing to K(x,x + L) come from Cl{~^^ + V) and 
Kd{L, L)Sl{-£i + V). By ([EM]), ([ESI), and 

K{x, x + L)^-l V(y) + i (^i ^,(2/) dy) 

which, given the periodicity of Vq, is (|5.15p . □ 
There is a second proof of (|5.13p . It depends on this analog of Naiman's lemma: 
Lemma 5.8. IfV is Lf^^ ^^^jj and + V commutes with Sl + S-l, then 

V{x + L) = V{x) (5.50) 
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Proof. Suppose first that V is bounded. Then 5*^ + S^l leaves D{—-^) invariant 
and commutes with it, so Sl + S-l commutes with V. If / is supported in a 
small neighborhood of a;o, (xq — (5, + 5) with \5\ < L/2, then {Sl + S-l){V f) is 
two separate pieces V(x — L)f{x — L) supported near xq + L and V{x + L)f(x + 
L) supported near xq — L, while V{Sl + S-l)/ is two pieces V{x)f{x ~ L) and 
V{x)f{x + L). Since the pieces are disjoint, 

V{x)f{x^L)^V{x + L)f{x-L) 

which implies (|5.50p . 

For general V, take g G C(j"(K) with J g{x) dx — 1 and note that / g{x)Sx{—-j^ + 

V)S-^x dx is —-^+9*y and it commutes with Sl + S^l also. But g^V is bounded, 
so it is periodic. (j5.50p follows by using an approximate (5-function. □ 

6. Block Jacobi Matrices and Matrix Orthogonal Polynomials 

What the magic formula suggests is that the Jacobi matrix J has parameters 
that approach an isospectral torus if and only if A( J) approaches + S^p. A( J) 
is a matrix of width 2p+l (i.e., A{J)m = if ^ {0, ±1, . ..,±p}) and Sp + S-p 
is a matrix with I's at the extremes. 

A matrix of width 2£ + 1 has the structure of a tridiagonal matrix if rewritten 
in terms of ^ x ^ blocks and + corresponds to Bk = 0, Ak = 1, the identity 
matrix, so A(J) ~ 5^ + S^^, at the matrix level, approaches the 'free case.' This 
will allow us to reduce our main theorems to matrix analogs of the theorems on 
perturbations of the free case. 

Of course, the association of the block matrix to orthogonal polynomials is 
critical — the orthogonality will be with respect to a matrix-valued measure. There 
is a huge literature on MOPRL (see, e.g., [MIMISIISIISIISSIIMIISZIISHIIMI 
[511 m [751 [m [119]) and MOPUC (see, e.g., [i[25l[2i[2l[2a[29l[30l[4l[70l[73l 
(M l I113| I120j ). In this section, our main purpose is to set notation and discuss the 
important notion of equivalent families of block Jacobi matrices, a notion discussed 
more explicitly in [21j . 

Given a semi-infinite complex matrix M = {mij}i<cij<oa and £ = 1,2,..., we 
define the £ x £ block decomposition as the family oi £ x £ matrices {Mqr}i<q r<oo 

by 

iMqr)tj = mi(^q-l)+i,l(r-l)+j (6.1) 

Definition. A block Jacobi matrix is an M where 



Mgr = 



Bq 


if r 


= g> 1 


Aq 


if r 


= q+l,q>l 




if r 


= q-l,q>2 





k- 


r\>2 



(6.2) 



with each Aq invertible and each Bq Hermitian. Here, following [102) Section 2.13], 
we use ^ for Hermitian adjoint; this is to avoid confusion with the Szego dual $* 
appearing in OPUC. 

We will start writing for such matrices. 

In analogy with the scalar case, one may be tempted to require 



Aq>0 



(6.3) 
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but to include A (J), we do not want to do that exclusively. If (|6.3p holds, we say 
that J is of type 1. 
If instead 

. . . v4„ > (6.4) 

for all ri, we say J is of type 2. 

An i'Kl matrix, if, is said to be in L, if it is lower triangular with strictly positive 
diagonal elements, that is, 

fo ifi<7 , , 

K,, = <^ ■' 6.5 

^ [ > if i = j ' 

If each Aq e £, we say that J is of type 3. 
The calculations in Section [3] and |4] show: 

Proposition 6.1. (i) //A is the discriminant oj a 'periodic Jacobi matrix, Jq, of 
period £, then for any Jacobi matrix, J, A( J) = is a block Jacobi matrix of 
type 3. 

(ii) // A is the discriminant of a periodic CMV matrix, Co, of even period £, then 
for any CMV matrix, C, of A{C) = J is a block Jacobi matrix of type 3. 

We will see that distinct j/'s may correspond to the same measure. Indeed, in 
the scalar case, the 6„'s and |a„|'s are fixed by the measure, but the arg(a„)'s are 
arbitrary. Thus, we define 

Definition. Two block Jacobi matrices, J and J , are called equivalent if and only 
if there is an £ x block diagonal unitary U = 1 J72 ffi C/3 © • • • (we will use Ui for 
1) so that 

J ^ UJU-^ (6.6) 

This is equivalent to 

B„ = C/„S„C/„-i (6.7a) 

i„ = \J^A,,U-l^ (6.7b) 

We will be interested in equivalence classes of J^'s. 

Proposition 6.2 ( 21J). Each equivalence class of J's has exactly one element 
each of type 1, type 2, and type 3. 

Definition. The Nevai class is the set of J^'s for which 

S„ AiA^ ^ 1 (6.8) 

The following is immediate from (j6.7p : 
Proposition 6.3. If some J is in the Nevai class, so are all equivalent J's. 

Damanik, Pushnitski, and Simon [2T prove that 

Proposition 6.4 ( 21 ). If J is in the Nevai class and is type 1, type 2, or type 3, 
then 

An ^ 1 (6.9) 
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We will sometimes need the MOPRL, the matrix orthogonal polynomials. What 
we describe here are the left OPs. There are also right OPs (see [21]), which we do 
not need here. An ^-dimensional matrix-valued measure is a positive scalar measure 
drjt{x) and a nonnegative I y. t matrix- valued function M{x). The matrix- valued 
measure 

dri{x) = M(x)drit(x) (6.10) 

can always be normalized by 

Tr(M(a;)) = t (6.11) 
We will always assume drj is normalized, that is, 

dT^{x) = 1 (6.12) 

The proper notion of nontriviality is a little subtle; it is discussed in detail in [21j . 
For our purpose here, it is sufhcient that ((•, defined below is nondegenerate on 
polynomials. 

If f,g are two ^-dimensional matrix- valued functions, we define 

((/,5))l = j g{x)M{x)f{x)Ur^t{x) 

Note that this 'inner product' returns matrix values. Recall also that ^ denotes the 
Hermitian conjugate of a matrix. The subscript 'L' is for 'left' and reflects the fact 
that if C is an X matrix, then 

((/,Cg))L = C((/,5))L (6.13) 
= ((/,. g))LCt (6.14) 

We will normally just write ((•,•)) from now on. 
Left orthonormal polynomials are of the form 

Pn{x) — Knx" + lowcr Order 

with matrix coefficients, defined by 

{(Pn,Pm)) = Snm^ (6.15) 

So long as drj is nontrivial, the p„ exist. They are not unique since if {J/njJ^^ 
unitary £ x £ matrices, 

Pn{x) = Un+lPn{x) (6.16) 

are also MOPRL. We demand kq = 1, that is, ^0(2;) — 1, and so [/i = 1. 

{pj}"=Q are a left module basis for matrix polynomials of degree n, that is, if / 
is any polynomial of degree n, then there are unique £ x £ matrices /□,...,/„ so 
that 

n 

f{x) ^^fjPjix) 

Indeed, 

fj^iiPjJ)) (6.17) 

For n = 1, 2, . . ., define 

Bn ^ {{Pn-l,xp„-i)) An ^ {{pn,XPn-l)) (6.18) 

Then, since xpj = X^fco '-'^P^ imphes {{pj,xpn)) = {{xpj,Pn)) = if j < n - 2, we 
have 

XPn{x) = An+lPn+l{x) + Bn+lPn {x) + Aip.n-i{x) (6.19) 
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This implies An+iKn+i — k„ so 

K„ = (Ai...A„)^i (6.20) 

and the type 2 condition is equivalent to k„ > 0. 

Looking at (|6.19|) . we see that (|6.16p holds for Pn,Pn if and only if An,Bn are 
related to An, Bn by (|6.7p . Jacobi matrix equivalence is just a 'change of phase' in 
the MOPRL. 

Given a block Jacobi matrix, we can view it as acting on the Hilbert space 
^2({1, 2, . . . }, C^) with inner product 

OO 

(/,.9) =E</"'5n)c^ (6.21) 

n=l 

If {ejYj^^ is the standard basis of C^, then {Sk-j}'kLij=i: defined by 

iSk;j)n = SknBj (6.22) 

is a basis. J acts on 2, . . . }, C^) via 

iJf)n = + Bnfn + A„/„+i (6.23) 

(with Aq^O). 

The spectral measure for J' is the £ x £ matrix-valued measure with 

{So;jJ{J)So;k) = J fix)dri,k{x) (6.24) 

for any scalar-valued function /. It is easy to see (e.g., [H]) that this map from 
J to rj inverts the one given by forming the MOPRL and defining J by (|6.18p . 
Moreover, and are equivalent if and only if dfj = drj. 
The m-function is defined by 

m{E) = J dvix) (6.25) 

^{So-.,iJ~E)-'So-.) (6.26) 

It is an £ X £ matrix-valued Herglotz function: 

Im£; > ImTO(£;) > (6.27) 

that is, -^{iTi — m^) is positive definite in the upper half-plane. For information on 
matrix Herglotz functions, see[2ll[Tll[Ta|4l|4i|4S[53[5l[Ml[M[63[9i HISl . 

Obviously, by ((05)) . m is constant over equivalence classes. 

As in the scalar case (see |102[ Section 1.2]), one has that for a.e. a; € M, 

lim^j^o m^x + ie) = m{x + iO) exists and 

dilac{x) = TT^^ Imm(a; + iO) dx (6.28) 

Here 

d77ac(2;) = M{x) d7]t;M (6.29) 

where dr]t-a,c is the a.c. part of drjt. Alternatively, drjuc is the unique matrix- valued 
measure which is a.c. (i.e., riac{I) = for any set with |/| = 0) and where there is 
a set K with \K\ = so (77 - ri^)(R \K) = 0. 

Given a block Jacobi matrix, J, by we mean the matrix with the top n 
(block matrix) rows and leftmost n columns removed, that is, 

= Sfc+„ Ai") = A,+„ (6.30) 
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We write m(")(z) for the m-function associated to J''^"-'. Equivalent J^'s do not have 
the same m'") for n > 1 (although they are unitarily related). We see m^^^ = m. 

We will need the following result of Aptekarev-Nikishin 'G' (see also \TD), a 
matrix analog of the well-known Jacobi-Stieltjes recursion for OPRL: 

Theorem 6.5 ([51 [H])- We have that 

m^^\E)-^ =E- B„+i - A„+im("+i)(£;)Ai+i (6.31) 
/orn = 0,1,2,.... 

Next, we need to note the following analog of a well-known scalar result (see, 
e.g., |104j ) proven in [2T] : 

Theorem 6.6. Let he a block Jacobi matrix with cross(i7) C [a, b]. Then, for any 
e, there is a K so that for k > K , 

cj{J^^'^)<z[a~e,b + e] (6.32) 

Finally, we need to look at poles and zeros of det(m(z)). In the scalar — 1) 
case, poles occur precisely at eigenvalues of J and zeros at eigenvalues of J^'"^', the 
once stripped Jacobi matrix. In that scalar case, these eigenvalues are distinct. 

In the matrix case, S and can have eigenvalues in common (as can be easily 
arranged by taking a direct sum of suitable scalar J's) so there can be cancellations. 
We say a scalar meromorphic function, f{z), has a zero/pole of order fc G Z at 
zo € C if (0 — zq)~^ f{z) has a finite nonzero limit as z — > zq. We will need the 
following result from [21]: 

Theorem 6.7 (,2lj). Let xq G M. Let qq be the multiplicity of xq as an eigenvalue 
of J , and qi its multiplicity as an eigenvalue of J\ . Then 

(a) go + gi < ^ 

(b) det(m(z)) has a zero/pole of order q\ — go- 

We will also need the following result from Aptekarev-Nikishin [^: 

Theorem 6.8. Let J be a block Jacobi matrix with (7oss(i/) — [—2,2], cr(i7) \ 
Ccssli/) a finite set and with g{x) — drji^dx) / dx we have 

y (4 - x^y^'"^ log(det(5(a;)) dx > -00 (6.33) 

Suppose J is type 2. Then 

lim Ai . . . An 

n — *oo 

exists and is a strictly positive matrix. 

7. A Denisov-Rakhmanov Theorem for MOPRL 

As preparation for proving Theorem ll.2l in Section[51 in this section we will prove 

Theorem 7.1. Let drj be a nondegenerate i x £ matrix-valued measure on R with 
associated block Jacobi matrix J of type 3 so that 

(i) aessC^) = [-2,2] (7.1) 

(ii) dr] = f{x) dx + drj^ (7.2) 

with drjs singular and 

det(/(x)) > (7.3) 
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-2,2]. Then 

Bn ^ An ^ 1 (7.4) 



Remark. ()7.3p says the a.c. spectrum has multipUcity i. 

If (17. 1|) is replaced by the stronger cr( J) = [—2,2] and type 3 by type 2, this is a 
theorem of Yakhlef-Marcellan [118j . We wiU prove Theorem 1 7. II by modifying their 
proof. 

The shift from type 2 to 3 is easy on accomit of Proposition 16.41 By applying 



the argument of |118j . we get yl„ — 1 for the equivalent J of type 2, conclude the 
whole equivalence class is in the Nevai class, and see An ^1. So we will only worry 
about the changes needed to go from a{J) = [—2, 2] to acss{J) — [—2, 2], where we 
follow Denisov's approach for the scalar case [31]. 

|118j relies on a matrix version of Rakhmanov's theorem proven by van Assche 
[113j . We need to extend it slightly to allow a.c. spectrum on a large subset of i9D 
rather than all of 9D: 

Theorem 7.2. Let dfj, be an £ x £ matrix-valued measure on and let {a„}^Q 
denote its matrix Verblunsky coefficients . Suppose 

d9 

d^i = w{e) — + dyi^ (7.5) 
2w 

where d^s is singular, and let 

n = {6: Aet{w{e)) > 0} (7.6) 

Then 

limsup |la„|l < 2V2^ (^1 - (^^^ ^ (7.7) 

Remarks. 1. For notation on MOPUC, see [21]. 

2. Where we use {an}5^Loi van Assche |113j uses {Hn}^=i related to a„ by 

Hn = -aLi (7-8) 

3. We follow notation from [113] and the variant of the scalar proof as in [103[ 
Section 9.1] where a„,5„,c„,(i„ below all appear. 

We define 



Cn,, = ^£^Tr(^^(^..)^L^^(^..)-l[^L^^(^..)t]-l^L(^..)t)l/2rf^ 

dn ^ ^J^\ri[^llie^')wie)^l:{e^y]'^')M 

Proposition 7.3. For every n > 0, we have that 

dn < bn,q for cvcry q>l (7.9) 

bl^g < 8f{l - Cn,q) for every q > 1 (7.10) 

dl < inf cn,g (7.11) 

9>1 
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Moreover, we have that 

/IOI\ 

<liminfd„ (7.12) 



2n 



Consequently, 



|f^|X3x 1/2 



limsup a„ < 2\/2£( 1 - ( -LJ- ) ) (7.13) 



Proof. The second to last displayed formula on |1131 p. 7] is (|7.9p . The estimates 
on the bottom half of |1131 p. 12] show that 

<'4 |l([^^(e^'')^J;+,(e^'^)-i][^^(e''')^J;+,(e'^)-i]t)V2_j|,2rf0 







21 

< — {ini - ine Cn q) 
= 8^2(1 - c„,,) 

which is ([77TO| . The third displayed formula on [TT3l p. 14] is (fTTTj) . 
Now, mimicking the estimates on the bottom half of [113|, p. 14], 

Tr([/(0)«;(0)/(0)t]i/4)d0 

< (^2^^^Tr(/(0)^L(^i.)-i(^L(^,.)t)-i/(0)t)^y^'(2^^rf„)i/2 
Taking n ~* oo, we see that 



< [2tt(Mv / /(6l)d/i(6l)/(6')t)i/4(27r£ liminf d„)^^^ 
V Jo. 

Removing the singular part as in [113j . we obtain 

i:T{[f{e)w{e)f{en'')d9 

n 

1/4 

1/2 



< 



27r£Tr / f{0)w{0)f{9)Ue\ (27r£ liminf d„ 



Proceeding as in [1131 pp. 15-16], it then follows that 

\n\£ < {2TTe\n\lY/'^{2TTt liminf d„) 

n — >OQ 



1/2 



which imphes (I7.12p . 

Putting these estimates together, 

a„ < bn,i < 2V2e{l - c„,i)i/2 < 2V2£(1 - ^2)1/2 

and hence 

limsup On < 2V2£(^1 - ^ 
which is dZHni). □ 
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In particular, for 27r — \Q\ small, 



lflW'\ 1/2 



limsup a„ = 0(^(^1 - (^yj j j =0((27r- 1^1)1/2) 

as in the scalar case. 

We have thus proven Theorem 17.21 To get Theorem 17.11 we follow |118j using 
the analog of Denisov's arguments for the case i = 1. 



Proof of Theorem \7.1\ By Proposition l6.41 we need only prove for the type 2 choice, 
for any e > 0, we have 

limsup(||i„-l|| + ||B„||) <e (7.14) 

By the Szego mapping and Geronimus connection formulae in [118j . this holds by 
Theorem 17.21 so long as for any e' > 0, we can find k so cf{J'^''^) C [—2 — e',2 + e'], 
and this is true by Theorem 16.61 □ 



8. A Denisov-Rakhmanov Theorem for Periodic OPRL 



Our main goal in this section is to prove Theorem 11.21 We will also prove the 
'hard' half of Theorem 11.11 The simplicity of the proof shows the magic in the 
magic formula! 



Proof of Theorem M.SX By a right limit of J, we mean a two-sided Jacobi matrix, 
Jr, (but with some a's allowed to vanish) so that for some subsequence rij oo 
and any k d Z, 

flnj+fc {0'r)k bnj+k ^ {br)k (8.1) 

By our standing convention, Jacobi parameters are uniformly bounded, so by com- 
pactness, if dm{{a, b),Tjg) 0, there exists a right limit Jr ^ Tj^. Thus, it suffices 
to show that any right limit Jr has Jr G Tjo • 

By the hypotheses of Theorem 11.21 the spectral mapping theorem, and the fact 
that A maps cress('/o) to [—2, 2] with a p-fold cover on (—2, 2), we see that 

A(J)ess = [-2,2] 

and A( J) has a.c. spectrum of multiplicity p. So thinking of = A( J) as a block 
Jacobi matrix, J' is of type 3 and the hypotheses of Theorem 17.11 apply. It follows 
that An 1, Bn 0. This means that 

A( J,) = SP + S"" 

so by the magic formula (Theorem 13. ip . Jr G Tj^. □ 

Rakhmanov's theorem is often related to issues of w-lim|p„p d/i and to the den- 
sity of zeros. We note that there are also results of that genre here: 

Theorem 8.1. If Jo is a periodic Jacobi matrix of period p and J is a Jacobi matrix 
with bounded Jacobi parameters whose right limits all lie in Tj^ [in particular, if 
the hypotheses of Theorem \l.S\ hold), then {with dji the measure for J) 
(a) 



w-lim - y \pjj^n{x)\ dfi{x) ~ diy (8.2) 

the density of zeros for Jq . 
(b) The density of zeros o/p„(x) converges to dv. 
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Proof. If Ji G Tjf, , then the spectral measure d^\^^^ associated to 5k has period p 

in k since Ji is periodic. Thus hniAr^oo 2W+1 H\j\<N ^A'j"''^ = | Si=i ^'^i"''^' 
the hmit is d;/ by the discussion in Subsection 12.81 Since 

lx%f\x)fdf,ix)^{S„.fS,) (8.3) 

and Jjj, Jj,j±i is very close to some {Ji)j.j, {Ji)j,j+i for |j — j'ol < M for fixed M 
and jo 00, we see that moments of LHS of (|8.2p are close to moments of dv. 
This proves (a). 

If J(") denotes the top left n X n submatrix of J, then 

/ x^d;.„ = iTr((j("))0 

so 




and thus (a) implies (b). □ 

We also have 
Theorem 8.2. // d„((a, 6), Tj„) 0, then 

(J) C (Jcss{Jo) 

Proof. By the magic formula, compactness, and the fact that every right limit of 
J is in Tjg, we see that every right limit of A(J) is + S~p, that is, An 1, 
B„ 0. Thus, by Weyl's theorem, (7ess(A(J)) = [-2,2]. Since 

(A(J))-A( 

we see a-css(J) C A"i([-2,2]) = (Tcss(Jo)- □ 
Remarks. 1. Since A is p to 1, we cannot conclude that acss{Jo) — o'css(«/) from 

A(CTcss(Jo)) = A( 

(J))- 

2. That crcss(^o) C (Tess(^) is a simple trial function argument given that J must 
have some right limits; see [551 [55] . 



9. Denisov-Rakhmanov Sets 

In this section, we want to show how one can take suitable limits of Theorem ll.2l 
to get a 'cheap' proof of similar theorems in other nonperiodic cases. We will 
also present an insight into the proper general form of Denisov-Rakhmanov-type 
theorems. 

The right limits we have discussed so far involve the weak product topology on 
the Jacobi parameters, so we will emphasize this fact by using the phrase 'weak 
right limits' in this section. We are also interested in limits in the £°°-topology for 
two-sided sequences, that is, {cn'^^}5^_oo ~^ {cn°°^}5^L_oo this topology if and 
only if, as fc — > 00, 

sup|cW-c(r)|^0 

n 

In terms of weak limits, we note the following: 
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Proposition 9.1. Let £ he a closed set. Let J be a half-line Jacobi matrix with 

Sac(J) =CToss(J) =f (9.1) 

Let Jr be a weak right limit of J . Then 

Sac( Jr) = <j{Jr) - £ (9.2) 

Remark. Note that (|9.2|) has a{Jr), not merely acss{Jr)- 
Proof. By resuhs in [68] . 

a{Jr) a £ C Sac(>/r) 

Since Sac(Jr) C a{Jr) trivially, ((9?2)) holds. □ 

Recall that a sequence {cra}5^_oo is called uniformly almost periodic (in the 
general theory of almost periodic functions, this defines 'almost periodic' — we add 
'imiformly' because the term is sometimes used in a weaker sense in the spectral 
theory hterature) if and only if {c'^^}^_3^ given by (c(^^)„ = Cn+e has compact 
closure in the ^°°-topology. 

Definition. A set £ is called essentially perfect if and only if £ is closed, and for 
al\E e£ and 6 > 0, \{E - d, E + 6) £\ > 0. 

Remark. Essentially perfect sets are precisely the sets, £, for which there is a purely 
a.c. measure dri with supp(c?r;) — £. 

Definition. A set £ is said to be a Denisov-Rakhmanov set if and only if 

(i) £ is essentially perfect and bounded. 

(ii) There is a set compact in the uniform topology so that for any bounded 
Jacobi matrix, J, for which (|9.ip holds, the set of right limits of J lies in Tg. 

The definition says nothing explicit about Tg being a torus, but by Propo- 
sition 19.11 if Jr e Tg, then ()9.2p holds, and since is closed under transla- 
tions, each Jr in is almost periodic. By Kotani theory (see [5^ [Ml I103j ). 
{Sn, {Jr — E ~ ie)^^5n) has real boundary values for a.e. E. In many cases and, 
in particular, if 5 is a finite union of closed intervals, Sodin-Yuditskii [108) (see 
also [5l I14j) proved there is a natural torus so that any almost periodic Jr with 
real boundary values lies in this torus. Thus for such cases, that 5 is a Denisov- 
Rakhmanov set can be connected to approach to an isospectral torus. In particular, 
our Theorem 11.21 implies the statement that (Jcaa{Ja) is a Denisov-Rakhmanov set. 

Given an essentially perfect set, £, we define 'D{£) to be the set of Jacobi matrices 
obeying (|9.1|) . 

The following two simple results will be the basis of our approximation theorems: 

Proposition 9.2. Let £ be an essentially perfect set. Suppose there are uniformly 
compact sets and T'^°°^ of two-sided Jacobi matrices so that 

(1) // J„ e T(") and J„ -> Joo weakly, then J^ e T^°°\ 

(2) For any weak right limit point Jr of some J G 'D{£), there is J & T*^"^ so 

sup |aW-a,H-|feW-6,|<i (9.3) 

\j\<n ^ 

Then £ is a Denisov-Rakhmanov set. 

Proof. Let J„ be the J guaranteed by (|9.3p . Then clearly, J„ converges weakly to Jr 
so, by (1), Jr G T(°°^. Since T(°°' is uniformly compact, £ is a Denisov-Rakhmanov 
set. □ 
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Proposition 9.3. Let Jq be a fixed periodic Jacobi matrix with essential spectrum 
So- Then for any n, there is a S > so that for any set £ with 

(a) £ <z{E: diiit{E, £o) < 6} (9.4) 

(b) \£\>{l-S)\£o\ (9.5) 

and any J e T^{£), we have that any right limit, Jr, obeys (|9.3p for some J G Tj^. 

Moreover, if p is fixed and C is a compact subset of [(0, oo) x W]p , then 6 can be 
picked to work for all Jq = {(a„, bn)}n=i G C . 

Proof. The uniformity claimed in the last statement comes from noting that choices 
can be made uniformly in the proof below. 

Let p be the period of Jq. We first claim that given Si, we can find S so if £ 
obeys (|01l-(l931). then 



dist(A(£), [-2,2]) < ^1 (9.6) 
\{x e (-2, 2) : aU p solutions of A{E) = x lie m£}\> A- Si (9.7) 

This is immediate from the continuity of A and its derivatives. 

Next, we note that given Si, we can find Si so that if J' is a p x p block Jacobi 
matrix so that 

'TessC^) C [-2^6^,2 + Si] 

\{E e [—2, 2] : J' has a.c. spectrum at E of multiplicity p}\ > 4 — Si 

then 

limsup\Jkm^{SP + S-P)km\<ei 

k,m — ^oo 

The proof of this is identical to the proof of the matrix Denisov-Rakhmanov theo- 
rem. 

Combining these steps, we are reduced to showing for any n and e, there is Si 
so for all two-sided J,- with dist(cr( J,.), £) < e, we have that 

sup \A{Jr) - {SP + S~P)k„^\ < ei (9.8) 

k,rn 

implies there is a J e Tjo so that (|9.3p holds. To do this, we first follow the proof 
of Theorem 13. II to note that for n, £2, and £3 fixed, we can find ei so (|9.8p implies 
there is a p-periodic such that 

||A(J«)-A(J)|| <e2 (9.9) 

and 

sup |af)-a«H-|6f)-6,|<£3 

|j|<" 

Finally, a compactness argument shows that for any n, we can find £4 so for any 
periodic j" with 

||A(J»)-(5^' + 5-f)|| <e4 

there is a J €E Tj^ so that 

ii^'-J|i<^ 

Putting these together implies l|9.3p . □ 
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Theorem 9.4. Let £±,£2, ■ ■ ■ be an arbitrary sequence in (2, 3, 4, ... ). For any £1- 
periodic Jacobi matrix J^^\ there exist A:2, fca, . . . so that for any limit periodic J 
with Jacobi coefficients 

00 

= a^°^ + ReiA^e^""/^!^-^"-] (9.10) 

00 

bn - fol"^ + Re[B„e2""/^--^"] (9.11) 

m=2 

obeying 

\A„,\ + \B„,\ < fc„ (9.12) 
we have that cr[J) is a Denisov-Rakhmanov set. 

Remark. The study of limit periodic discrete Schrodinger operators with smah 
tails was initiated by Avron-Simon [7j and Chulaevsky 19]. They prove purely a.c. 
spectrum. 

Proof. As in [HHn], one can pick the fc^'s so the spectrum is purely a.c. and so 
that the union of all isospectral tori for the periodic approximates lie in a fixed £°° 
compact sets. This implies the limit periodic potentials also have compact isospec- 
tral sets, and within this compact set, weak convergence implies norm convergence 
so hypothesis (1) of Proposition 19.21 holds. By decreasing the fcm's if necessary, 
Proposition 19.31 continuity of the spectrum in norm, and absolute continuity 
of periodic spectrum imply we can be sure that (|9.3p holds. Thus Proposition 19.21 
implies this theorem. □ 

Our final theorem in this section is the following: 

Theorem 9.5. Fix £. Let G = {(ai, 02, . . . , a^+i, /J^+i) G M.'^^+'^ : ai < Pi < 
a2 < /?2 < ■ • ■ < f3i+i}. For {a, (3) € G, define 

i+i 

£(dj)^[j[a„l3,] (9.13) 

Then {(a, /3) : £{d,j3) is a Denisov-Rahkmanov set\ contains a dense Gs- 

Remarks. 1. As we have seen, the £{d,P) which arise from periodic problems are 
precisely those where the harmonic measure of each Cj = [aj,(3j] is rational. In 
particular, if we fix a and /?£+!, the set of /3's that are periodic is countable, and 
so certainly not a Gs- We show that the family that leads to Denisov-Rahkmanov 
sets is uncountable. 

2. It is a reasonable conjecture that every £ is a Denisov-Rakhmanov set, so 
this result is weak. We include it because it is such a 'cheap' way to go beyond the 
periodic case using only that case. 

Proof. For each {a, (3) , it is known [108] that there is an isospectral torus T of 
almost periodic J's where (whole line) spectrum is precisely £{d,(3). It follows 
from the construction in [105] that if (a("),/3(")) G G converge to {d^°^\ (3'-°°'^) £ G, 
then condition (1) of Proposition 19.21 holds . 

Let Gp be the subset of G coming from periodic problems — this is dense in G- 
For (a("),^(")) e Gp, pick J{d'^°\ 0'^°^) periodic with f (a(o)^^(o)) as spectrum and 
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pick 6n{d^°\ P^^'') via Proposition 19.31 requiring (5„ < ^min(|/3]"' — aj°^|,|a^'^j — 
/3f^|). Let t/(")(ao,A)) = : E{dJ) obeys (EU/dnSl) for £ = £{dj), 

So = £{a^°\f]^°'^) and 6 = (5„}, and let 

C/(") =U;7(")(a(°),/3(°)) 

This is dense and open. Then n„C7(") is a dense whose points, by construction 
and Proposition 19.21 correspond to Denisov-Rakhmanov sets. □ 

10. Sum Rules for MOPRL 

In this section, our main goal is to prove the following two theorems about block 
Jacobi matrices: 

Theorem 10.1 (P2 Sum Rule for MOPRL). Let J' be a block Jacobi matrix with 
i X £ Jacobi parameters {yl„}5^]^, {_B„}^]^ and matrix measure drj. Let m{E) be 
given by (I6.25P and suppose cross(i7) ~ [~2, 2]. Define for z(^B)\{z = E + E^^ : 
i?ea(J)\[-2,2]} 

M{z) ^ ~m{z + z-^) (10.1) 

Let F, G be the functions 

F{(3 + r')^W-r^~\og(3^] (10.2) 
/or /3 G M \ [-1, 1], that is, = /3 + /J-^ e M \ [-2, 2] and 

G{a) ^ ~l~\og{a^) aG(0,oo) (10.3) 
Then Van^^i M{re'^^) exists for a.e. 9 and 



h /^"g( det(irJ(e-^)) ) 

00 

+ F{E)^Y.^,{\Bl + \G{\A^\)) 

Ee<j(J)\[-2,2\ n=l 



(10.4) 



Remarks. 1. All terms are positive (since F and G are positive, this is evident for 
two terms; positivity of the integral will be seen below), so this sum rule always 
makes sense, although some terms may be +00. 



2. Recall that |A„| = y AJiA„; although since the formula for G{a) only involves 
a^, one does not need to take a square-root. 

3. Because of the trace and absolute value, Tr(-|i?^ + -iGd^nj)) is constant over 
equivalence classes of Jacobi matrix parameters. 

4. In the type 1 case, the RHS of (fTa4| is finite if and only ii J - Sp - S-p is 
Hilbert-Schmidt. This is also true when J' is of type 3; see Proposition 1 1 1 . 1 2l 

Theorem 10.2 (Sharp Case Co Sum Rule for MOPRL). Consider the three quan- 
tities: 

1,1 sm 



Z{J) = fog -r-r. — irrr-jp:^ dO (10.5) 

^ ' At: Jq ^Vdet(ImAf(e^''))y ^ ' 

SoiJ)^ J2 log(l/3|) (10.6) 

E<faiJ)\[-2^2] 
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where j3 is related to E by 

^eM\hl,l] E^f3 + p-^ (10.7) 

and 



N 

AoiJ)= lim -V log(det(|A„|)) (10.8) 

N^OQ ^ — ^ 

n=l 

which we suppose exists but it may be +oo or —oo. Then 

(i) // any two of Z,£q, Aq are finite, then so is the third. 

(ii) // all are finite, then 

Z{J)^AoiJ)+£o{J) (10.9) 

(iii) // all are finite, then 

N 

lim y Tr(B„) (10.10) 



N~ „_ 

n=l 



exists. 



Remark. We will prove (and actually use it to prove Theorem 1 1.3p that if £q{J) < 
oo, then Z{J) < oo so long as 

AaiJ) = liniinf ('-y log(det(|A„|))') < oo (10.11) 

Theorem llO.ll is a matrix- valued analog of the OPRL P2 sum rule of Killip-Simon 
[5T] . and Theorem[Tn2]of the OPRL Case Co sum rule by Simon-Zlatos tlOTj. Both 
were refinements of sum rules of Case |17[ 118] who in turn was motivated by earlier 
KdV and Toda sum rules. Case only considered short-range |a„ — 1| -I- |6„|, while 
[5T1 1107j considered the necessary techniques to go up to the borderline of validity. 
[107] had some simplifications of [BT], and [100| further simplified, although each 
of the later two proofs depends heavily on the earlier ones. Here, following Simon 
[100] . we will prove a nonlocal step-by-step sum rule. As there, the key is a suitable 
representation theorem for meromorphic Herglotz functions — in this case, extended 
to matrix- valued functions. 

For a S (—1,1), we define Blaschke factors as usual by 

{°-^ < a < 1 
-l<a<0 ^ ' 

l — az — 

Proposition 10.3. Let f{z) be an £ x £ matrix-valued meromorphic function on D 
so that 

(i) ± Im/(z) > w/ien ±Iniz > (10.13) 

(ii) lim f{z)z-'^ = 1 (10.14) 

where liR f = ^{f — f^) . Then 

(a) For a.e. 9, limrti /(re*") = /(e'") exists. 

(b) log|det(/(e*"))| Gni<p<„oin<9©,d^/2^) 

(c) All the zeros and poles of det(/(z)) lie on (—1, 1) and are of finite order. 
Let {zjyj^i and {pj}°°^i be those zeros and poles of det{f{z)) repeated up to 
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multiplicity {it can also happen that both sets are finite), z = is not included 
in {zj}. Then 

Boo{z) = hm (10.15) 

exists and obeys: 

(i) Boo j-s analytic and nonvanishing on C\{zj}Ll{pj}Ll{zJ^}Ll{pJ^}Ll{±l} 

(ii) |Boo(e*«)| = 1 ondB\{±l} 
(iii) 

|arg(Boo(z))| < 2Tr£ (10.16) 

for \z\ < 1 with arg normalized by argi?oo(0) = 0. 
(d) We have the representation 

det(/(z)) = z'B^{z)cxp(^l log|det(/(e»^))| ^) (10.17) 

Remarks. 1. It should be possible to prove that < arg(i?oo(^)) < ti"^ for Imz > 0; 
we settle for the weaker result. 

2. (|10.14p is not central for a result of this type, but it is true in applications 
and simplifies the notation. 

3. This result for ^ = 1 is in [lOOj . i > 1 has some subtleties, but the basic 
strategy we use is that of [lOOj . 

We will prove this result in a sequence of lemmas: 

Lemma 10.4. det(/(z)) is analytic and nonvanishing m = {z : z € D, Imz > 
0}, and arg(det(/(z))) can be chosen in that region to be continuous so that 

< arg(det(/(z))) < 7r£ (10.18) 

Proof. In Q, all matrix elements {ip, f(z)ip) are analytic and have a.e. boundary 
values (since they are scalar Herglotz functions), so by polarization, f{z) is analytic 
on rt and has a.e. boundary values. Thus det(/(z)) as a polynomial in matrix 
elements is analytic on n. 
Consider 

P(A,z) = det(Al-/(z)) (10.19) 

which is a polynomial in A with analytic coefficients away from the poles of /. It 
follows, for z near any zq about which / is analytic, that the roots P(A, z) ~ 
written as a function of z are analytic functions in (z — zg)^/*^ for some k depend- 
ing on Zq. It then follows that near any fixed zq, all roots are analytic, that is, 
singularities are isolated. 

Pick xo £ (0,e) so that xof{xQ) > 0, so all eigenvalues Ai(a;o), . . . , A£(a;o) are 
in (0, oo). Let z e 17 be a point about which all eigenvalues are analytic, and 
let 7(z) be a simple closed path from xo to z which avoids the discrete set where 
eigenvalues are not analytic and lie in f2 except for xq with, say, 7(0) = xq, 7(1) = z. 
By analytically continuing eigenvalues, we get function {Aj (z)}^^]^, so Aj(z) are all 
the eigenvalues of f{j{t)) and Xj{0) € (0,oo). By Im/ > 0, ImAj(z) > 0, so if we 
define arg(Aj(z)) with arg(Aj(0)) = 0, we have 

< arg(Aj(z)) < tt 



42 



D. DAMANIK, R. KILLIP, AND B. SIMON 



Thus 

£ 

arg(det(/(z))) = ^arg(A,(^)) 

normalized by arg(det(/(xo))) = obeys (|10.18|) . 

By analyticity of det(/(z)) and the fact that it is nonvanishing, arg(det(/(z))) is 
uniquely defined as a continuous function on f2 with limej^o a'rg(det(/(a;o + — 0. 
By the above, (jlO.lSp holds at all points ^ in $7 where all eigenvalues are analytic 
and so, by continuity and the open mapping theorem for analytic functions, all 
points. □ 

Lemma 10.5. Let a < b lie in (—1, 1) so that both a and b are neither a zero nor 
a pole o/det(/(z)). Let Z(a,b), P{a,b) be the number of zeros, poles o/det(/(z)) 
in (a, b) counting multiplicity. Then 

\Z{a,b)-Pia,b)\<£ (10.20) 

Proof. By the argument principle, 27r(Z — P) is the change of arg(det(/(2;)) along 
the circle through a and b centered at ^ (a + 6) . By Lemma 110. 4[ this is at most 
2{en). □ 

Lemma 10.6. The sets of zeros and poles {with multiplicity) o/det(/(z)), including 

(k) (k) 

the £-fold zero at z ^ 0, can be written as £ subsets z^ ,Pj with k = and 
—Nk < j < Nk {with Nk and Nk among 1,2,..., oo) so that Zq'"'' — and 

zf < pf^ < zf^, (10.21) 

for all allowed values of j . 



Remarks. 1. If there are infinitely many z in (—1, 0) and in (0, 1), then Nk = Nk = 
oo for all k. The awkwardness requiring Nk, Nk is only needed if there are finitely 
many zeros. 

2. To avoid notational complexity, we slightly lied if Nk or Nk is finite. If Nk is 

(k) (k) 

finite, Zj runs to j = Nk- Pj can then run to either Nk or Nk — 1. 

Proof. Construct Si, S2, ■ . ■ , Si as follows: Set Zg^"* — 0. Let pg^'' be the first pole 
larger than Zq^\ z[^'^ the first zero larger than p'^\ p^p the next pole, etc. This 
either continues indefinitely, in which case we set A^i —00, or stops because there 
is no next zero or pole. Then do the same to the left of 0, that is, p*^j is the first 

pole smaller than z^ , etc. Clearly, the points in 5*1 obey (|10.2ip . Now remove 
the points of 5*1 (or decrease their multiplicity by 1) and repeat the construction 
(starting with Zg — 0) to make S2, S^, . . . , St. 

We claim that after we construct £ Sj^s, we have exhausted all the poles and 
zeros. Let us show this is true for (0,1); the argument for (—1,0) is similar (and 
since has multiplicity £, it is removed after £ steps). 

Suppose z is a zero that is left and it is closer to zero than any leftover zero 
or pole. If z lies in some (p^'^^ zj^\), j — 0, 1, ... , we could have used it as zj^\ 
so it cannot lie in any such interval. Put differently, there are only matched zeros 
and poles in (0, z) n Uj^^Sj. By the choice of z, there are no other poles in (0, z). 
Thus, for small S, the interval {—6, z + S) has £+1 extra zeros over poles, violating 
Lemma [10.51 So the closest leftover point is not a zero. 
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Suppose p is a pole that is left and it is closer to zero than any other leftover 
zero or pole. As above, p cannot lie in any (zj'^^p^'^''), j = 0, 1, . . . , so there are 
only matched zeros and poles in [0,p) n Uj^iSj. But then, for small S, {S,p + S) 
has £+1 extra poles, violating Lemma [10.51 Thus U^j^^Sj includes all zeros and 
poles. □ 

Lemma 10.7. The limit Boo{z) of pO.lSp exists and obeys conditions (i)~(iii) of 
Proposition \1 0.3( c) . 

Proof. Renumber the p'j'^ into a single sequence pi,p2, • ■ • , so \pi\ < \p2\ < • • • 
and let Zm be the correspondingly paired ^j+i (paired to the p^*^'' that is p) . Since 
{ipf \ are disjoint subsets of (0, 1) for each fixed k, 

oo oo 



3 = i 



so we see that 



Y,\^,-p,\<2i 

The existence of Sqo then follows by Proposition 13.8.2 of [103 , as do (i) and (ii). 

To get (iii), we note that just taking the zeros and poles in a single Sj yields a 
set obeying (13.8.5) and (13.8.6) of [103j . So, by (13.8.10), that product has arg 
bounded by 27r. The £-fold product thus obeys (|10.16p . □ 

Proof of Provosition \10.S[ Given Lemma [10.71 the proof is essentially that of The- 
orem 13.8.3 of [TUg] . 

, , _ det(/(z)) 
'^'^ = ^^W^iz) 

is analytic and nonvanishing on B) with g{{)) > (since Boo{Q) > 0). Moreover, by 

pim)i and (nnni), 

|argg(z)| < 47r£ (10.23) 

so, by M. Riesz' theorem, \og{g{z)) G np<oo^^^(9D) from which (a), (b) of the 
theorem are immediate and (d) follows from the Poisson representation for \og{g{z)) 
since log(|5(e'^)|) = log(|det(/(e'«))|. □ 

Now we turn to block Jacobi matrices where we obtain: 

Theorem 10.8 (Nonlocal Step-by-Step Sum Rule for Block Jacobi Matrices). 

Let be a block Jacobi matrix with (Toss(i/) C [—2,2] and Jacobi parameters 
{AmBn\^=i. Let J''^^^ denote this Jacobi matrix with the top row of blocks and 
left-most column of blocks removed. Let m{E) , m^^\E) be the m-functions given by 
(j^:^ . Let M, M(i) be defined on D by 

M{z) ^ -m{z + z-^) (10.24) 

with poles at {pi}fLi where pi + p^^ are eigenvalues of J . We repeat each pi a 
number of times equal to the multiplicity of the eigenvalues {equivalently, the rank 
of the residue) . Let{zi}^J[^ be the corresponding points for J^^^\ Then 
(a) The Blaschke product, Boo{z), defined by the {zi\yj {pi\ via ()10.15p exists and 
obeys (i)-(iii) of Proposition ] 10. "M c). 
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(b) M{z) and Af(i)(2) have limits M{e''^) and M'^^'>{e''^) as r 1 for z = re'" for 
a.e. 9 in and 

det(ImM)(e'^) 



log 



det(ImA/(i)(e*»)) 



l<p<oo 



2^ 



det 



\Ai\M{z] 



Boc{z) exp 



e'" + z 



log 



det(ImM(e*'')) 



det(ImA//(i)(e*«)) 



d9 
An 



(10.25) 



(10.26) 



Remark. As in the case i = 1, it can happen (although not in examples where 
sum rules are finite) that det(ImM(e*^)) = det(Im M(i)(e**)) = for 61 in a set of 
positive measure, (b) and (c) are shorthand for the more precise 

(i) For a.e. 9, det(ImM(e*^)) = if and only if det(ImM(i)(e*^)) = 0. 

(ii) There is an a.e. positive function (7(6*) on9D, equal to det(ImM(e*^))/ det(Im Af*^^)(e 
when the ratio is not 0/0 so that (|10.25p and (|10.26p hold if the formal ratio 

is replaced by g{9). 

Proof. Given Proposition 110. 3i this is essentially identical to the proof of Theo- 
rem 13.8.4 of [103] with care given to matrix issues. We begin by noting that (|6.3ip 
for n — !■ n + 1 first implies near z = 



and then by ([QS)) that 

'M(")(z)^ -1 



M("+i)(z)-i = + 0(1) (10.27) 



= 1 - Bn+iz - (4+i^n+i - 1)2' + 0{zJ) (10.28) 



Since M(")(z)/z is near 1 for z small, we can compute its determinant using 



det(C) = exp(Tr(log(C)) 
which holds if \\C - 1|| < 1. Thus 

log det ( 

= Tr(i3„+i)^ + [Tr{([4_,iA„+i ~ 1] + ^BI^,)}]z' + 0{z^) 
In addition, (|6.3ip implies 

Im[M(z)"^] = Im(z + z"^) - Ai 1\tvMi{z)A\ 
so at points where M{z) has radial limits (a.e. 6*, see below), 

- [M(e''')^]-ilmM(e''')[M(e*^)]-i = -Ai Im Af(e*^)4 
which, using (on account of dct(|C|)2 ^ det(C^) det(C)) 

|det(Ai)| =det(|Ai|) 

yields 

|, WM \Af( ^o^^\2 dct(ImM(e'^)) 
|det(|^i|M(e ))| =det(IniA/i(e^«)) 

We now apply Proposition ll0.3l to M{z) which obeys (|10.13p (since Im(z+z^^) < 
on ID) and (|10.24p has a minus sign) and (|10.14p by (|10.27p . 



(10.29) 

(10.30) 

(10.31) 
(10.32) 



(10.33) 
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By Theorem 16 .71 our Boc{z) here (after perhaps cancehng some zeros and poles) 
is the Boo{z) of Proposition 110.31 (a) and (b) immediately follow from Proposi- 
tion [TH^l We get (fTU:^ from (fTUT7|) by using (noting (fTU:^ has a \/Att 
while (|10.17p a l/27r on account of the square on the left side of (|10.33p ). We also 
use that if c is a positive constant, 

-pf/^log(c^)?Uc (10.34) 



z Att 



□ 



As in |100j . we can get step- by-step P2 (originally in [6T]), Co, Ci (originally in 
[107] ) sum rules immediately from Taylor expansion of the log of (|10.26p . We let 
fiji<7) be the numbers in (—1, 1) \ {0} for which Ej = (3j + are eigenvalues of 
J' counting multiplicities. 

Theorem 10.9 (Co, Ci, P2 Step-by-Step Sum Rules), 
(i) 

J_ / det(ImAf(i)(e'^)) Y 
AttJ, det(ImAf(e''^)) J 

-^log(|/?,(^)|)-log(|/3,(:7«)|) = -log(det|Ai|) 



(ii) 



1 r^'^, /det(ImM(i)(e'^)) 
log' 



27: Jo det(ImAf(e»»)) , 

+ T.i^AJ) iPi'iJ)-')] mj^'^) ~ mj'^'''))-'] ^ Tr(Si) 

(iii) 

1 /det(ImAf(l)(e»»))\ . ^^n^^^a^ 

AttJ \ det{ImMieJ'>)) j ^ '^"^ ' (10.37) 
where F is given by (110. 2p and G by (jl0.3p . 

Remark. The Ej{J) in (-00, -2) and (2, 00) and Ej{J^^^) interlace in the ^ = 1 
case. In the general £ case, we have at most I fewer eigenvalues of J^^^^ on any 
(—00, —Eq) or {Eq, 00) so, as in Lemma |10.6[ one can decompose into £ interlacing 
subsets. This and the monotonicity of functions like F show the eigenvalue sums 
in (I10.35p - (|10.37p are conditionally convergent. Similarly, the integrals are always 
convergent. 

Proof. Apply log to both sides of p0.26p and take Taylor coefficients. The constant 
term is (|10.35l) and the first derivative is (|10.36p . If L{z) is the log of the left side 
and R{z) of the right, then 

L{0) + lL"{0) = R{0) + ^R"{0) 

is ([10371. □ 
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The proofs of Theorems I10.1H10.21 are now identical to those of the scalar case; 
see, for example, the discussion of Theorems 13.8.6 and 13.8.8 of |103j . In partic- 
ular, Z{J) and Q{J) (the integral on the right of (|10.4p ) are negatives of relative 
entropies, and so, lower semi-continuous. 

11. SZEGO AND KiLLIP-SlMON THEOREMS WHEN AlL GAPS ArE OpEN 

Our goal here is to prove Theorems 11.31 and 11.41 Our strategy, of course, will be 
to translate Theorems 110.11 and 110.21 for A(J) to statements about J. Firstly, we 
need to relate the a.c. part of the matrix measure for A(J) to the a.c. part of the 
(scalar) measure for J. And secondly, to relate P norms of coefficients of A (J) to 
the distance of J's Jacobi parameters to the isospectral torus. We begin with the 
first question. Thus we take 

dT]j{x) = uj{x) dx + d7]j,s{x) (11-1) 

with dfij^s singular and uj supported precisely on o'css(^o)- By this assumption and 
the spectral mapping theorem, A( J) has a.c. spectrum precisely on [—2, 2] so the 
matrix measure for A (J) has the form 

dVAiJ)iE) = W{E)dE + d7]^^j)^,{E) (11.2) 

Proposition 11.1. Let Jq be a periodic Jacobi matrix with period p and J a Jacobi 
matrix with Jacobi parameters {a„, fonj^^i and measure drjj of the form (|11.1|) with 
LO supported on acsa{J())- Let A he the discriminant for Jq and W{E) the a.c. part 
of the p X p matrix-valued measure dri^(^j-^ associated to A(J) (so W is a p x p 
matrix). Then for E G (—2,2) and A^^{E) — {xi, . . . ,Xp}, 

detiw{E)) = (nap) Yfi«r)Tn^(^^)) (i^-^) 

Proof. In the block Jacobi form, drj^(^j) has jk matrix element equal to the spectral 
measure of the operator A(J) associated to Sj,Sk, that is, j F{x){dri\{x))jk = 
{Sj,F{A{J))Sk). But Sj =pj_i(J)^i. It follows that 

p 

WkjiE) = ^t^(x,)(|A'(a:,)|)"V-i(a;£)p,-i(a;^) (11-4) 

£=1 

Note that the factors of 1/A' arise from the Jacobian ^ = A'{x). We can re- write 
pi.4p as Wkj{E) = {MAM*)kj where A is the diagonal matrix 

Aim = SeraC^{xe){\A'{xi)\y'^ (11.5) 

and M is the matrix 

Mke=Pk-i{xi) k = l,...,p; i=l,...,p (11.6) 
Next we compute det(M); dot (A) is easy. Note that 

.k-l s^l 

Pk-i{xi) = I ]^ I -I- lower order (11-7) 

Moreover, inductively one sees that the lower order terms can be neglected in the 
determinant — they can be removed by subtracting a multiple of rows above (i.e.. 
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smaller values of k). Thus, 

r p / fe-i \ -In 

det(Af) = 



1.-1 V i-i / J / S>fc 



by the well-known formula for Vandermonde determinants. This can be simplified 
further. The points Xj are precisely the zeros of the polynomial A(a;) — E; hence, 
invoking (|2.14p . 

Aix)-E=(flafA 'ff[(^-x,)\ 
In this way we discover that 

det(Af)2 = (n«r) Yn«f )Tn i^'(^^-)i) (i^-^) 

Multiplying this by det(A) gives (ITO)) . □ 
Corollary 11.2. //Jo /las all gaps open and a > — 1, then 

J (A- E^)°'\logdet{W {E))\dE < oo (11.10) 



i/ and only if 

dist{x,R\(T^ss{Jo)r\^oguj{x)\dx < oo (11-11) 



siJo) 

When a — — ^, the same conclusion holds even if some gaps are closed. 

Remark. Since a > —1, (4 — E'^)" (resp., dist(. . . )") are in for some p > 1, so 
the log_,_ ( ) is always integrable and these conditions are equivalent to the integral 
without I 'I being larger than — cx). 

Proof. Changing variables via E = A(x) and applying Proposition 1 1 1 . ll shows that 
(|11.10|1 holds if and only if 

I logw(x)| (4 - A(a;)2)° \A'{x)\dx < oo (11-12) 

= (Jo) 

If all gaps are open then |A'(a;)| is strictly positive on acss{Jo), while 4 — A(a;)^ 
is a polynomial with a simple zero at each band edge (and no others). This proves 
the first claim. 

At a closed gap, 4 — A(x)^ has a double zero and A'(x) a simple zero. When 
a — these cancel exactly. □ 

Next we turn to the issue. Given any two-sided periodic matrix J with Jacobi 
parameters {a„,fen}^=i and fixed periodic Jq, let Bjg{J), Ajg{J) be the constant 
pxp blocks in Ajg{J). We are heading towards showing that ||-Bjo(J)||2 + II^Jq(^)~ 
III2 is comparable to dist((a„, fe„)^^;^, T/q). This will be the key to showing £^ tails 
in the matrix pieces of Ajg(J) — — S^^ for general J is equivalent to (|1.28|) . changes 

The crucial fact will be that the polynomial coefficients of Ajg{J) — Ajq(Jo) 
are comparable to dist((a„, bn)n=ii'^Jo)- ^'^^ this we need the following, which is a changes 
simple application of the implicit function theorem and compactness: 
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Lemma 11.3. Let F be a C°° map of an open set U C M" to with £ < n. 
Suppose T = F^^{yo) is a smooth manifold of dimension n — £ and compact for 
some yo M.^ , and 

rank((VF)(a;o)) = £ (11.13) 

for all xq G T. Then for any compact neighborhood, K, of T , there are CK,dK G 
(0, oo) so for all x G K , 

CK\F{x)-yn\ <dist{x,T) <dK\F{x)~yn\ (11.14) 

One can restate (jll.l3p in a more illuminated way in terms of the components 
Fi, . . . , Fg oiT. Of course, VF^ (xq) is orthogonal to T at xq. The condition (|11.13|) 
is equivalent to saying that {W jF{xo)}j^i span the normal bundle to T. This is 
equivalent to saying they are linearly independent. Notice that if Jq has all gaps 
open, Tjg is of dimension p — 1 — 2p — (p -\- 1) and Aj^, is a polynomial of degree p, 
hence with p + 1 coefficients. Thus the following shows we can use Lemma lll.3l 

Theorem 11.4. Suppose all gaps are open for some periodic Jq. Then at any point 
in T/,-, , the gradients of the derivatives of the coefficients of Aj span the normal 
bundle of Tj„ m M^p. 

Proof. Aj(, has the form 

p p 
Ajo(a;) = (ai . . . 0^)"^ W{x - Xj) = ^ CjX^ 
i=i 3=0 
where \j are the roots. The coefficients thus obey 

ciCp^ ^ 51 Afei . . . Afei_j = Sfe_f £<p 

l</ci<---<fcp_£<p 

It is well known that if 

p 

ti^Y.^] (11.17) 

then tg is £si plus a polynomial in {sj}^zi[, so {Vijjj^]^ and {Vsj}^^]^ span the 
same space. It follows that we need only show the gradients of c" ^ and if span the 
normal bundle of Tj^ . 
Let 

Mo - {(a„, 6„)Li : aia2 . . . flp = af\f^ . . . af^ ■ b^ + ■ ■ ■ + bp ^ fe^"^ + • • • + 

(11.18) 

We know Ai^ D Tjg. Clearly, Vc^ ^ and Vti span the normal bundle to Mq since 
ti = J2^=i ^3 (see (|2.14p ). Thus we need only show the projections of {Vt^j^^j 
into the tangent space of A^o span the normal bundle of Tj^ in A^o- 

Studies of the Toda flows show that A4o is a symplectic manifold with {ti}g^2 
Poisson commuting. Since the symplectic form on Mq is nondegenerate, to say 
{Vij}^^2 span the normal bundle is the same as saying that the Hamiltonian 
flows generated by {tj}^^2 span the tangent bundle of T/p, or equivalently, given 
dimiTjg) — p — 1, that these Hamiltonian flows are independent. 

This independence is a theorem of van Moerbeke [114[ Theorem 5.2] or [106] . □ 



(11.15) 
(11.16) 
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Lemma 11.5. Let Xk be the projection onto the k- dimensional space spanned by 
{^j}j=i- For any compact subset, K, of period p Jacobi matrices, there exist con- 
stants ck and dx in (0, oo) so for all J G K , 



CK 



< d 



K 



(11.19) 



Proof. {J^Xp+i}^=o s-re independent since has strictly positive elements in the 
£-th diagonal and {J''}k<i only has zero elements there. Hence, the matrix 

Tr(xp+iJ^j'^Xp+i)|,,^o_ (11.20) 

is strictly positive so (|11.19p holds for each fixed J. The optimal constants are 
clearly continuous so uniformly bounded above and below on K . □ 



Proposition 11.6. Let Jq be a periodic Jacobi matrix with all gaps open. For any 
compact neighborhood K of Tj^ in (0, oo)p x W, there are constants ck and dx in 
(0, oo) so that for all J Cz K , 



ck{\\Aj,{J) 



Proof. We have that 



+ \\Bj,{J)f)'^'<dist{J,Tj„) 

<dK{\\Aj,{J) 



IIP 



\BjoiJ)\ 



2U/2 



2\\Aj,{J) - If + \\Bj,{J)r < ||[A,;„(J) - (5^ + S-n]Xp\\ 



But by the magic formula, 



<4Pjo(J)-] 



so 



Xp+i 



where Ci is the difference of coefficients for J and Jq . By Lemma lll.5| 



\\[Aj,ij)-isp+s-n]xp+ir-Y. 



(11.21) 

(11.22) 
(11.23) 

(11.24) 



£=0 



where ~ means the ratio is bounded above and away from zero on compact subsets. 
By Lemma [11. 31 and Theorem lll.4[ 

p 

^|c,|2^dist(J,Tjj2 (11.25) 

Combining this with pi.24p proves the proposition. □ 

Now we take a general J not periodic and form Ajp(J) which is a one-sided 
block Jacobi matrix with block elements Anjg{J),Bnj„{J). 



Lemma 11.7. Ajq(J)m for k < £ depends only on {bjjjZk^a '^"'^ {O'jYjtk-i 
where a = \_-^{p — {£ ^ k))\ is the greatest integer less than or equal to ^[p—{i-'k)]. 

Proof. Each factor of J changes index by at most one. In order to get from k to £, 
£ — k steps are needed. The remainder cannot go below ^ — a or above k + a and 
get back to fc in p steps. □ 
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Lemma 11.8. Let J have Jacobi parameters {fln, 6n}j5°=i • Liet J be periodic with 
period p and suppose &„ = 6„ for kp ~ p < n < kp + 2p and a„ — a„ for hp ^ p < 
n < kp + 2p ^ 1. Then 



Ak,Jo (J) = Ajo (J) Bk,j„ (J) = Bj„ (J) 
Proof. Immediate from Lemma 111.71 



□ 



Lemma 11.9. Let k < £ and a = [i(p — — k))]. For any two J and J and any 

K , there is Ck so that 



\^Jo{J)m- ^Jo{J)kA<CK sup [\bj~b,\ + \a,-~aj\] (11.26) 

k — a.<j<.l-\'a. 



SO long as 



^M\bj\ + \bA + \a,\ + \a,\]<K 



(11.27) 



Proof. Immediate from Lemma 111.71 and the fact that Aj^ has matrix elements 
that are fixed (given Jq) polynomials in a's and 6's. □ 

Lemma 11.10. (a) For any Jacobi matrix, J , and £ — 1,2, . . m = 1,2, . . 

(11.28) 

and for £ = 2,3, . . ., m = 1,2, . . ., 



H.— L 



j=0 



(b) For Jq periodic of period p >2 and m — 1,2, . 

' ' ' ^m+p— 1 



(0) (0) 



(11.29) 



(11.30) 



A,/o(J)mm+p_l = [a'fj ■ ■ -anj+p^i) (flm ' ' ' am+p-2) (^^(^m+j " b, 

3=0 



do) - 

in+j. 



(11.31) 



Proof, (a) Since J changes index by at most one, 

{^J )m m-\-£ — {Jmm+l) ■ • • {Jm+e-lm+i) 

proving (|11.28p . while 



J=0 



m+j m+£— 1 



which, given (lll.28p . proves (|11.29p . 
(b) By 



p-i 

3=0 



p-2^ 



A,„(J) = (4°)...a("))-^ 

which, given (a), {,P'~^)mm+p = {J^^^)mm+p-i = if A: = 2,3,..., and the 
periodicity of a'^^^ and yields (|11.30p and (|11.31|) . □ 
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Lemma 11.11. Suppose that Ajg{J) — — S ^ is a Hilhert- Schmidt operator on 
^2({0,1,2,...}). Then, 

^(a„a„+i • • • On+p-i - ai^^a^l^^i • • • ai°jp_i)^ < oo (11.32) 



5^ 5^(fe„+,-6i"i.) <oo (11.33) 

n \j=Q j 

^(a„+p - a„)^ < oo (11.34) 

n 

^(&„+p-fe„)' <oo (11.35) 

n 

Proof. For a Hilbcrt-Sclimidt operator, any subset of matrix elements lies in so 
by (fTT^ni , 

^ |a„ • • • a„+p_i(a^°) • • • a^°l^_^) ^ - 1 < oo 

n 

which, given that al"'' • • ■a''^^p_^ is n- independent, impUes (|11.32p . 

Similarly, (|11.3ip implies (jll.33p if we note that {aj} bounded and a„ ■ ■ ■ fln+p-i — 
af'"' ■ ■ ■ a^p"^ > implies inf aj > 0, so 

inf (a^^ • • • a|°L_i) ^(a^ • • • am+p-2) > 

Since the difference of sequences is i'^ , (111.32^ implies (since an^ is periodic) 

f<oo 

n 

which, given that inf Oj > 0, implies (|11.34p . 
Similarly, since 

p-i 

^^(^n+l+j ~ ^n+j) = bn+p — 6n 
3=0 

(fTOS]) imphes (fTL35l) . □ 

Our next preliminary is to relate A ^ C to 

\A\ = yCdU (11.36) 

Proposition 11.12. The map A i— > from C to positive definite matrices is a 
diffeomorphism. In particular, for A's in C with — 1|| < i, there exist constants 
Ci and C2 so that 

Ci\\A- III2 < II |A| - III2 < C2IIA- III2 (11.37) 

Proof. By (|11.36p . A 1-^ |A| is a smooth map. 

The inverse map (strictly \A\'^ ^ A) is known as the Cholesky factorization; 
see [52l I115| . Given i? > 0, apply the Gram-Schmidt procedure to the (linearly 
independent) columns of B working from right to left. This gives a factorization 
B = QA with Q unitary and A ^ C. Note that |yl| — B and that because the 
columns of B are linearly dependent, B i-^ A is also a smooth map. □ 
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Theorem 11.13. Let Jq be a two-sided p-periodic Jacobi matrix with all gaps open 
and let Aj^ denote its discriminant. For a Jacobi matrix with parameters (a„,6„), 
the following are equivalent: 



(i) Ajo(J) - SP- S^P IS a Hilbert-Schmidt operator on ^^({o, 1, 2, . . .}). 

(ii) E„Tr{i32 + |A„-l|2}<oo. 

(iii) E„Tr{B2 + (|A„|-l)2}<oo. 

(iv) E„Tr{B2+G(|A„|)}<oo. 

(v) Y..mdm{{a,b)Xlo? 

(vi) Y.md.^{{a,b).Xlo? 



Here we have adopted the abbreviations An '■— Anj„{J) and Bn '■= Bn.Jo{J) 

Proof. amounts to the definition of the Hilbert-Schmidt norm. 

(ii) <;4>(iii) follows from Proposition 1 11.1 21 

(iii) <;4>(iv) Notice that G, defined in p0.3p . obeys 

ce{x - 1)^ < G{x) <c',{x-lf V X e (e, e~^) 

Applying this to the eigenvalues of |A„| yields this equivalence. 

(v) <^(vi) is the q = 2 case of Proposition 13. 51 

(vi) =»(i) By Lemma ril.7[ each matrix element of Aj^ (J) — Sp — S~p is a smooth 
function of p consecutive pairs (a„,6„); moreover, by the magic formula, all of 
these smooth functions vanish if the corresponding p-tuple belongs to Tjg. The 
implication now follows from the fact that smooth functions are Lipschitz. 

(i)^(vi) Define J^'"'-' to be the p-periodic Jacobi matrix that equals J on block 



for ^ = 1, 2, . . . Obviously, jC^) = J on block k and, by pTM)) and (fTT35)) . the 
difference on blocks k — 1 and k + 1 are in £'^, that is. 



fc, that is. 




— Okp+e 



(11.38) 




(fc-l)p<j<(fe+2)p-l 



(11.39) 



Together with Lemma [11.91 and Proposition 1 11. 61 this implies 




(11.40) 



k 



On the other hand, (111.39P implies that for j — 1, . . . ,p, 




(11.41) 



k 



By the triangle inequality. 



dkp+j{J,Tj,f < 2dkp+j{J, J^^^f + 2dkp{J^''\Tj,f 



so (111.40P and (111.4ip imply (vi). □ 

Proof of Theorem \1.4\ We will refer to the three statements (i)-(iii) of Theorem 1 1.41 
simply by their numbers. 



Suppose first that 



d,„((a,6),rjj) (1142) 
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Then (i) holds by Theorem 11.11 Moreover, by Theorem 111.131 and the hypothesis 
that all gaps are open, the RHS of ()10.4|) is finite. Therefore, the LHS is finite. 
Next we use this fact to prove (ii) and (iii). 
As A' is nonvanishing at all band edges, 

N 

J^Fi^iE,)) < oo ^ ^dist(£;,,aess(J))^'" < oo (11.43) 

which verifies (iii). By Corollary 1 11.2[ 

Leftmost term in (|10.4p < oo -<=^ (n) holds (11.44) 

This completes the proof of (i)-(iii). 

Conversely, if (i)-(iii) hold, then by (|11.43p . (|11.44p . and (jTOl)) . we see that the 
RHS of pHi)) is finite. By Theorem [nH this implies (|11.42p . □ 

Proof of Theorem ] 1.3[ Let Pj be the /3's associated to A(J), that is, \(3j\ > 1, 
Pj+Pj^ = Ej with the eigenvalues of A( J) in M\ [-2, 2] . Then \og\(3j | ~ | - 1 
as P^±l smah and \(3j \ - 1 - {\Ej \ - 2)^/'^. Therefore, 

(Hnnxoo ^ ^(|£;,|-2)i/2 ^ (imD (11.45) 
j 

By Corollary [nH 

(fTHSll < 00 (fL^ (11.46) 

Finally, if {An, Bn}^^i are the p x p blocks in A(J), we have An — U\An\ for 
some U with |det(J7)| = 1, so 

n I. 

(11.47) 



(0) 



(0) 



deti^„i = idet(^„)i = n n 

j=(n-l)p+l k=j 

by (1331) and ^J^. Thus 

N Np 

5^1og(det(|A„|))-p^log 

n=l k=l 

is bounded. Thus (fL27l) is equivalent to (|10.11|) . 

By Theorem [1021 we see that when (fL24| holds, then (fL25| <^ (fr27|) . and if 
they hold, (|1.26p has a limit. 

Moreover, if they hold, the hypotheses of Theorem 11.41 hold, so (|1.28p is true. 
That (|1.29p holds is a theorem of Peherstorfer-Yuditskii [87l ; see also the remark 
below. □ 

In the remainder of this section, we will describe an alternate approach to proving 
(|1.29p ; one based on combining the magic formula with Theorem l6.8l Unfortunately, 
because of the strong hypothesis on the discrete spectrum that appears in this 
theorem, we will not recover the full formulation from Theorem 1 1.31 

Let denote the (unique) type 2 block Jacobi matrix that is equivalent to 
J = Ajq(J), which is of type 3. Further, let us use Aj and Aj to denote the 
off-diagonal block entries of J' and J, respectively. 

If we strengthen the hypothesis (|1.24p to finiteness of the discrete spectrum 
(i.e., finiteness of the set a{J') \ aess{>J)), then Theorem 16.81 shows that (|1.25p 
implies the convergence of the product Ai ■■■ An as n oo. In view of (|6.7bp and 
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Proposition 1 1 1 . 1 2l this convergence is inherited by the product Ai ■ ■ ■ An- Thus, it 
remains only to connect the convergence of this matrix product to the behavior of 
the sequences of parameters. This is the job of the next lemma. In the original 
version of this paper, it was only proved that the sequence {a„} was asymptotic 
to a fixed periodic sequence. The argument for the sequence {b„} was provided by 
one of the referees; we are most grateful for this. 

Lemma 11.14. Let Jq be a p-periodic two-sided Jacobi matrix and let A = Aj^ de- 
note its discriminant. Let J be a one-sided Jacobi matrix with parameters {a„, 6„}. 
Suppose the product A„ ■ ■ ■ Ai converges to a non-singular matrix as n —> oo . Here 
An and Bn denote the p x p block entries of A(J). Then the parameters of J 
asymptotically converge to fixed periodic parameters in the sense of (|1.29p . 

Proof. By applying the same affine transformation (i.e., x i— s- ax 4-/3) to both J and 
Jo, we may assume that the discriminant of Jq takes the form A{x) = x^ -\-0{xP~'^). 
This will significantly simplify some of the formulae that follow. Note also that this 
transformation makes a^^"^ ■ ■ ■ ap"-* ~ 1 and fe^"-* + • • • + b^'^ — 0. 
Let An and i?„ denote the p x p block entries of A (J). Then 

{An)k,k — ap(n-l)+fc ' ' ' Ipn+fc-l 
{An)k+l,k = ap(n-l)+fc+l ' ' ' Opn+fc-l [&p(n-l)+fc+l + ' ' ' + 6pTi+fc] 

as can be read off from Lemma lll.lOl Using this and the lower-triangular structure 
of the matrices Aj , one may quickly deduce 

pn+fc— 1 

{Ai---An)k^= n «J (11-48) 

j=k 



k,k 



{Ai ■ ■ ■ An)k+l,k — • • • {Ar-l)k+lM+l{Ar)k+l,k{Ar+l)k,k ' ' ' i^n) 

(11.49) 



r=l 
/pn-\-k—l 





To see that the sequence n i-^ Opn+k converges for each fixed fc e {1, . . . one 
need only take ratios of (|11.48p for consecutive values of k and the same n (and 
also for (n, k = p) and (n-\-l,k — 1)), then send n — > oo. 

For the parameters &„, one may proceed in a similar fashion: For example, when 
2 < fc < p — 1, the fact that 

{Ai - ■ ■ An)k+l,k {Ai ■ ■ ■ An)k.k~l , , 

{Al---An)k,k {Ai ■ ■ ■ An)k-l,k-l 

shows US that bpn+k converges as n — )■ oo. □ 



12. SZEGO AND KiLLIP-SlMON THEOREMS WhEN SOME GAPS ArE CLOSED 

Here we want to examine what might replace Theorems 11.31 and 11.41 if Jq is 
periodic but with some closed gaps. The Szego-type theorem is almost the same as 
Theorem O 
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Theorem 12.1. Let Jo be any two-sided periodic Jacobi matrix with Jacobi pa- 
rameters {an \ fen'''}^_oo7 '^^'^ J one-sided Jacobi matrix with Jacobi parameters 
{a„,6„}^]^ and spectral measure drj. Suppose that p.22p holds, and that 

N 

J2 dist(i;™, (Tess( J))'/' < OO (12.1) 
m— 1 

{Em}m=i is a labelling of the eigenvalues of J outside (Tcss{J)- Then 



log 

.{Jo) 



implies 



(^^^dist{x,R\a^ss{Jo)y^^^ dx < oo (12.2) 



pN , 



exists and lies in (— cxd,oo). Conversely, (jl2.2p holds so long as 

lim sup (^J2\og(^^^^ > -oo (12.4) 

and then the limit in (jl2.3p exists and lies in (—00,00). 

Moreover, if ()12.2|) or ()12.4p holds, then there is Ji S Tj^, so 

d™(J,Ji)->0 (12.5) 

Remark. All that is missing is (jl.28p which we do not claim. However, since 
p^ / P^ imply (i)-(iii) of Theorem [TU below, we have 

Proof. As noted, even with closed gaps, (|1.25p is equivalent to (jll.lOp (see Corol- 
lary [TT21)- Once one notes this, the proof of Theorem 11.31 provides all the results 
stated as Theorem II 2 .11 □ 



Theorem II .41 used open gaps in two ways. First, in the translation of a matrix 
pseudo-Szego condition, (|11.10p with a = i to the original spectral measure of 
J, and second, translating a Hilbert-Schmidt bound on A(J) — — to 
approach to the isospectral torus. The second issue can be finessed if we leave 
things as a Hilbert-Schmidt condition, which reduces to a sum of translates of 
an explicit positive polynomial in the a^'s and 6„'s being finite. As for translating 
()11.10p with a = ^ , the argument that proved CoroUary l 1 1 . 2l translates immediately 
to 

Lemma 12.2. Suppose (t(Jo) has closed gaps at {yj}j^i C cr(Jo). Then (|11.10p 
holds with a = ^ if and only if 



t 

/ dist{x,m\a{Jo)Y/^Y[\x-yj\^\log{uj{x))\dx < 



00 (12.6) 
Plugging this into our proof of Theorem 11.41 immediately yields 



Theorem 12.3. Let Jq be a two-sided periodic Jacobi matrix with closed gaps at 
{yj}^j^i C cr(Jo) and let J be a Jacobi matrix. Then 

Ty{{Aj,{J) - SP - S-Pf) < 00 (12.7) 



if and only if 
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(i) (fL22ll holds. 

(ii) (11.2411 holds with ^ replaced by 

(iii) IT^ holds. 



Example 12.4. Take Jo to be the two-sided free Jacobi matrix but think of it as 
period 2. Then 

Aj{x) ^x^ -2 

and a direct calculation of ,P — — S^^ shows that (|12.7p is equivalent to the 
three conditions 



< CO 



E 



E(""" 



n+1 



lY < oo 



(12.8) 
(12.9) 
(12.10) 



If bn = and 



(12.11) 



a„ = l + (-l)"(n+l)-'3 
then P^ - (|12.10p hold if and only if > i while, of course, ([11^ requires 
/? > i. This is one of many known extensions of the (Jp-free) Killip-Simon theorem 
(see, e.g., Laptev et al. 67 , Kupin 66J, and Nazarov et al. [81]). Some of these 
results have MOPRL analogs which, via the magic formula, lead to variants of 
Theorems O and [T231 



13. Eigenvalue Bounds for MOPRL 

There are Birman-Schwinger kernels for MOPRL and it should be possible to 
extend the proofs of most bounds on the number of eigenvalues outside [—2,2] or on 
moments of \Ej \ — 2 from the scalar to matrix case with optimal constants. But if 
one is willing to settle for less than optimal constants (but still not awful constants), 
there is a simple method to go from the scalar to matrix case. It depends on the 
following: 



Theorem 13.1. Let J be an £x£ block Jacobi matrix in the Nevai class with Jacobi 
parameters Let E^{J') denote its eigenvalues counting multiplicity 

outside [-2, 2], that is, > > ■ ■ ■ > 2 > -2 > 
scalar Jacobi matrix with a„ = 1 and 



> E^ > E^ 



±||B„||±M„-i-l||±||A 



1 



Let J± be the 
(13.1) 



and let J^£' be an £-fold direct sum of J± . Then 



\EfiJ)\ < \EfiJ^^ 



(13.2) 



Proof. The matrix analog of the observation of Hundertmark-Simon [55J extended 
to 2£ X 2i matrices (with £ x £ blocks) says that 



-\An-l\ 
1 



-\An-l 



< 





Ai 



An 





< 



I A. - 1| 



1 

\An-l 



since 



(0 c\ 


2 


fc^c 


^ 


[c^ o) 




\ 





|C|| 



(13.3) 



(13.4) 
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Thus writing Jj_ a,s £ x £ blocks with each block a multiple of 1, 

ji^^ <J< (13.5) 
from which p3.2|l is immediate. □ 
Corollary 13.2. For any block Jacobi matrix, J , in Nevai class, 

Y^iEfiJf <2£J2 \\Bn\\+MY, \\An-l\\ (13.6) 

j',ib n n 

Remark. In particular, this implies if the RHS of p3.6p is finite, so is the LHS. 
Proof . Hundcrtmark-Simon 55J proved 

^(i?±(J±)2- 4)1/2 <^5± (13 7) 

j ri 

from which (fT^TB)) follows by □ 

14. The Analog of Nevai's Conjecture 

Proof of Theorem[rM By (f02l) . J-Jq is trace class. Thus J'^-J^ = EI^o J'"{J- 
Jo)J^^^^'' is trace class, so Aj„ (J) - Aj„ ( Jq) = Aj„ (J) - {S^ + S^p) is trace class. 

It follows that if (|1.32p holds and A( J) has matrix Jacobi parameters {An, i?„}^i 
that 

oo oo 

I] l|S„|| +^ ||A„- 111 < oo (14.1) 

n— 1 n— 1 

By Corollary [1321 the eigenvalues A (J) obey 

oo 

^(|i?±|- 2)1/2 <oo (14 2) 

(|14.ip also implies 

oo 

^|log(det|^„|)| <(X3 (14.3) 
We can apply Theorem 110.21 and conclude that Z{J) is finite, that is, 

J {4- £;2)-i/2 log(det(Ty (£;))) dE > ~oo (14.4) 
By Corollary [TTH we obtain pT^ . □ 

15. Perturbations of Periodic OPUC 

In this final section, we want to present the translations of our results to OPUC. 
Since the magic formula maps periodic OPUC to MOPRL, the changes needed in 
the proofs will be minor, although for the analog of Theorem 11.41 there is one sig- 
nificant change. It is interesting to note the sequence of mappings for the OPUC 
periodic Rakhmanov's theorem. We map OPUC to MOPRL using the magic for- 
mula and them map that to MOPUC using the Szego map. 

The OPUC version of Theorem 11.11 is already in Last-Simon [68 . As for Theo- 
rem HH 
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Theorem 15.1. Let Cq be a two-sided periodic CM V matrix. Let C be an ordinary 
CMV matrix with Verblunsky coefficients {cknji^Lo- Suppose 

Sac(C) = cross(Co) (15-1) 

Then, as m —s- oo, 

dm(a,TcJ^O (15.2) 

Proof. Let Cr be a right limit of C. By TheoremO A{Cr) = Sp + S-P. Thus, by 
Theorem mH C,. e Tc,,. □ 

As for the analog of Theorem 11.31 if we drop discussion of i'^ convergence, it 
holds, similar to Theorem 112. II 

Theorem 15.2. Let Cq be a two-sided periodic CMV matrix with Verblunsky co- 
efficients {aj'^''}^_oo; and C a one-sided CMV matrix with Verblunsky coefficients 
{o^AfLo o.nd spectral measure dfj,. Suppose that CToss(C) = cr(Co) and 

N 

dist(£;„, ae,,(C))i/2 < oo (15.3) 

m— 1 

where {Ejn}^^! '^^ ^ labelling of the eigenvalues of C outside crcss(^o)- Then 

-j log(^)dist(0,R\aees(Co))-i/2^<oo (15.4) 

implies 

exists and lies in (—00,00). Conversely, (|1.25p holds so long as 



N 



and then the limit in (jl5.5p exists and lies in (—00,00). 

Remarks. 1. We have not stated that the q;„ have a limit in Tc^. We suspect that 
the methods of 87j extend to the OPUC but have not checked this and they do 
not explicitly mention it. 

2. Of course, if this theorem is applicable and Co obeys the conditions of Theo- 



rem 115.31 below, then we have a result on P convergence to . 
3. One can replace p"^^ by the logarithmic capacity of cr(Co). 

Proof. At open gaps, A'(e*^) 7^ 0, so (|15.3p is equivalent to 

^ (|i?|- 2)1/2 < 00 (^5 7) 

£^[-2,2] 
Be(T(A(C)) 

Moreover, by (|44l) -(l47 |) . we have that 

np j+P^^ 



log(det|A„|) log(det(yl„)) ^ log 



Pk_ 

j = (n-l)p+l k=j fk 



n n ^ 



(15,8) 



Now just follow the proof of Theorem 1 12. II □ 
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In carrying over Theorem 11.41 to OPUC, one runs into a serious roadblock: van 
Moerbeke's theorem |114) that the Hamiltonian flows generated by the coefficients 
of the tj (given by (|11.17p ) are independent is not known for OPUC. Instead, we 
use a weaker result of Simon |103[ Section 11.10] that proves the derivatives of 
coefficients of A span the normal bundle for a dense open set that is 'most' of the 
points with all gaps open. We will call the isospectral tori in this dense open set 
the generic independence tori. Then by mimicking the arguments in Section [TH we 
obtain 

Theorem 15.3. Let Co he a two-sided CMV matrix in a generic independence torus 
with Verblunsky coefficients {a^^^}'^_^^. Let C be a CMV matrix with Verblunsky 
coefficients {q;j}°^q. Then 

oo 

J2 dm{a,Tc„ f < ^ (15.9) 

if and only if 
(i) Cross(C) = (7 {Co) 

(ii) For the eigenvalues {Ej}jL^ not in a{Co), 

N 

^ dist(£;j, cr (Co ))^/^ < OO (15.10) 

(iii) // /i is the spectral measure for C, then 

'I logf%^Vist(^.5ID\aess(Co))i/'^ <oo (15.11) 
J<t(Co) \ d.0 J 27r 

Our last result is a periodic OPUC version of Nevai's conjecture. 

Theorem 15.4. Let Co be a two-sided p-periodic CMV matrix and letC be a CMV 
matrix with Verblunsky coefficients {o:j}°°^o- Then 

oo 

J2 dm{a,Tc„) < (X (15.12) 

m=0 

implies (115. 4p 

All the above results assume the period p is even. However, by using sieving (see 
Example 1.6.14 of ,102 ) to map periodp to period 2p, one can extend Theorems ll5.1l 
and 115.2] to p odd. In particular, we obtain the p = 1 results of pUl [5] as very 
special cases of ours. 
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